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ABSTRACT 


A high resolution, multi-level, primitive equation ocean model is used to examine 
the response to wind forcing of an idealized, flat-bottomed oceanic regime on a f-plane, 
along the eastern ocean boundary off the west coast of the Iberian Peninsula. Three 
experiments are conducted to investigate the role of different types of prescribed wind 
forcing. In experiment 1, a band of equatorward winds which are uniform alongshore 
but contain zonal variability (1.e., large negative wind stress curl near the coast), results 
in an equatorward coastal surface current nearshore and a poleward surface current 
offshore. With time, the currents become unstable and anticyclonic warm core eddies 
develop in the region of negative wind stress curl. In experiment 2, the model 1s forced 
with poleward wind stress for 20 days followed by a sudden change to equatorward wind 
stress. A 1-2 day transition from downwelling and a poleward surface current near the 
coast to upwelling and an equatorward surface jet occurs. With time the upwelled 
isotherms are displaced offshore, the equatorward surface jet widens and gets shallower 
nearshore, and the poleward undercurrent moves closer to the surface (depth ~ 120 m 
at the coast). For experiment 3, a time series of 6-hourly values of the north/south 
component of the wind, computed from Fleet Numerical Oceanography Center synoptic 
surface pressure analyses are used to investigate the response of the coastal current sys- 
tem to time dependent wind forcing. Equatorward wind events increase the vertical ex- 
tent and cross-section of the surface equatorward jet and reduce the poleward 
undercurrent. Relaxations reduce the vertical extent of the equatorward surface coastal 
current, increase the cross-sectional area of the poleward undercurrent, and the under- 
current tends to shoal nearshore. Poleward winds displace the equatorward current off- 
shore and create a poleward surface current nearshore. By the middle of the upwelling 
season, eddies are generated. By the end of the upwelling season (with the seasonal re- 
versal to poleward winds), the eddies and equatorward surface current weaken or disap- 
pear and a poleward surface current develops nearshore. The results from the 
experiments support the hypothesis that wind forcing, and particularly time dependent 
winds, can be a significant generation mechanism for eddies, jets, fronts and upwelling 


filaments off the west coast of the Iberian Peninsula. 
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I. INTRODUCTION 


A. OBJECTIVES 

Satellite images of the sea surface temperature distribution off Portugal and north- 
West Spain, taken during conditions favorable for upwelling, show several 
plumes/filaments extending offshore near 43° N, 41° 40’ N, 40° 20’ N, 38° 40’ N, and 
off southwest Portugal (Fiuza and Sousa, 1989). Ongoing studies, based on multi-year 
sequences of satellite imagery, indicate that such filaments tend to show up near the 
Same places, mainly between the submarine canyons which indent the continental mar- 
gin (Fiuza and Sousa, 1989). Observations have also shown that pairs of anticyclonic 
and cyclonic mesoscale eddies, with horizontal dimensions of order 100 km (Fiuza, 1984) 
can also exist off the west coast of Portugal and Spain. In other eastern boundary cur- 
rent regions of subtropical gyres, particularly in the California Current System, eddies 
have been shown to be generated by baroclinic instability (Emery and Mysak, 1980; 
Thomson, 1984; Ikeda and Emery, 1984; Batteen er a/., 1989), barotropic instability 
(Wright, 1980), or mixed barotropic/baroclinic instability (Thomson, 1984; Wright, 
1980) of the large scale flow. 

The primary objective of this thesis is to investigate the role(s) of climatological and 
time dependent wind forcing as possible generation mechanisni(s) for the coastal jets, 
undercurrents, fronts, and eddies in the eastern boundary current region off the west 
coast of the Iberian Peninsula (IP). Experiments, designed to reproduce some of the 
possible important and/or unique features of the study area, will be systematically run 
using a multi-level, eddy-resolving primitive equation (PE) ocean model developed by 
Haney (1985) and modified by Batteen er al. (1989). The water temperature field will 
be initialized in all experiments based on climatological data for the area. Different wind 
forcing will be used to study some of the important processes due to wind driven circu- 
lation for the area. The role of wind forcing in the generation of mesoscale features in 
the current system off the IP has not been systematically investigated and may be the 
most important mechanisin. Satellite infrared imagery has shown evidence of cold water 
filaments, jets, and eddies during periods of winds favorable for upwelling. These obser- 
vations provide evidence for wind forcing as a possible important mechanism for the 
formation of the mesoscale features. Eddies and jets could respond to the seasonal mean 


wind or to wind fluctuations with time scales of a few days during the upwelling season. 


The organization of the thesis is as follows. In Chapter I, the wind fields, hydrology, 
and circulation that are thought to be important or unique features for the study area 
are described. The PE ocean model and the experimental conditions used in the study 
are then presented (Chapter II). Chapter III describes the basis and procedures used to 
applv spectral analysis for the evaluation of wave lengths and kinetic energy of eddy 
fields that result from the simulations. In Chapter IV the results of the model simu- 
lations along with an analysis of the results are presented. A comparison of model re- 
sults with available observations is also made. Conclusions and recommendations for 


future work are presented in Chapter V. 


B. BACKGROUND 
In the following section, general descriptions of the wind fields, water masses, and 
circulation for the study region are presented. An introduction to the wind forcing used 
in some of the experiments 1s also given. 
1. The Wind Fields 

The seasonal regime of the Azores high pressure system and the meridional mi- 
gration of the subtropical front, which separates the mid-latitude from the tropical air 
masses, regulates the mean wind conditions off the IP west coast (Fiuza, 1982). The 
center of the Azores high migrates along 38° W, from 27° to 33° N between March and 
August. The mean east-west pressure contrast between the coast of Portugal and the 
high center is higher in summer (8 mb) than in winter (1 mb), resulting in weak mean 
westerly winds in winter, and considerably stronger northerly and northwesterly mean 
winds in summer (Fiuza, 1982). The northern limit of southward upwelling favorable 
winds shifts from about 27° N, near the Canary islands, in January to about 43° N, over 
the northern limit of the Portuguese west coast, in July (Wooster et al/., 1976). Based 
on one-degree latitude/longitude climatological summaries of wind stress estimates from 
ship reports near the coast of the IP, Bakun (1987) concluded that the mean wind stress 
curl distribution is weak and irregular during fall and early winter, but strongly cyclonic 
during the spring and summer, when the spatially separated wind stress curl maximum 
appears adjacent to the northwestern coast of the Iberian Peninsula and off southern 
Portugal (Figure 1). The region of cyclonic Wind stress curl extends from the coast up 
to about 200 km offshore; farther offshore there is anticyclonic wind stress curl (Bakun, 
fost): 

Local wind patterns near the coast, associated with orographic features such as 


capes and/or rias (fjord tvpe of estuary along the northwest coast of Spain), may signif- 
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Figure 1. Wind stress curl distribution for June-July in the area off the west coast 
of the IP and NW Africa: Units are dynes/cm’? per 1000 km. Regions 
of anticyclonic wind stress curl are shaded and regions of cyclonic wind 


Stress curl are unshaded (from Bakun, 1987). 


icantly influence the flow near the coast. For instance, McClain et al. (1986) combined 
information on the wind stress derived from surface pressure charts and ship observa- 
tions for the period from 18 to 22 April 1982 and found consistently strong winds with 
a heading to the southwest at the mouth of Ria de Arosa, parallel to the axis of the ria, 
and generally weaker winds a small distance offshore (about 20 km). The existence of 


the resulting large negative wind stress curl, possibly orographically induced, has been 


concluded by McClain er al. (1986) to be the primary driving mechanism for the current 
shear associated with an observed southward coastal jet adjacent to a northward off- 
shore current near the shelf break. These wind and current features will be explored in 
experiment |. 

A time series of monthly values (over a 32-year period, 1948-1979) of the wind 
stress in four, five-degree latitude/longitude quadrangles off the IP was analyzed by 
Bakun (1987). The quadrangles used in the analysis are shown in Figure 2. The area 
covered by these data overlays the numerical model domain used in the present study 
(see Chapter II). Based on his analysis of climatological and time-dependent winds, the 


following features of the wind field were noted by Bakun (1987): 


e The seasonal cycle (long term mean monthly values) of the alongshore wind stress 
in the near-coastal areas, shows strongest equatorward values during the upwelling 
season, roughly from April through September. 


e There is a general trend towards increased values of the alongshore (southward) 
wind stress (averaging values for April through September of each year) from the 
beginning of the time series (1948) to the end (1979). 


e The long term mean monthly values of alongshore wind stress are stronger in the 
southern areas compared to the northern areas. This feature 1s apparent if every 
month 1s included in the averaging process, and particularly so, if the upwelling 
season alone is considered. 

A time series of 6-hourly values of the north/south component of the wind ve- 
locity computed from Fleet Numerical Oceanography Center (F NOC) surface pressure 
analvses, for the period from | February to 31 December 1986, at 39°N, 10°W 1s shown 
in Figure 3. The wind velocity computation is based on a finite difference calculation 
of the geostrophic wind (from the analyzed surface pressure field interpolated onto a 
3-degree latitude/longitude grid), with a subsequent multiplication by 0.7 and a rotation 
of 15° to the left to account for friction. These wind values mav be seen as an approxi- 
mation to the large spatial scale, time dependent winds in the IP. Fluctuations of the 
alongshore component of the wind stress with time scales of a few days (which may be 
associated with reversals of the coastal currents) can be seen. In addition, sequences of 
alternating upwelling favorable events (southward winds) and relaxations (northward 
winds) are common features of the time series. Spectral densities computed from this 
data series (Figure 4) show a considerable amount of variance in the event time scale 
(e.g., 4-20 days). The strong peak at the diurnal frequency is probably associated with 


differential day/night heating over land and over the ocean and corresponding changes 





Figure 2. Locations of 5-degree latitude/longitude quadrangles off the west coast of 


the IP: (areas labeled ‘A’, ’B’, ‘C’, and ‘D’) used to summarize maritime 
reports in monthly time series. The mean position of reports within each 
of these areas, over the period 1948 to 1979, is marked by the ‘X’ svm- 
bols surrounded by small circles. The mean latitude and longitude of the 
reports from all four areas together (i.e., from the entire 10-degree 


quadrangle) is marked by the large asterisk (from Bakun, 1987). 


in the cross-shore pressure gradient. The response of currents to this time series will be 
explored in experiment 3 (see Chapter IV). 
2. The Water Masses 
Climatological temperature, salinity, and potential density data, from Levitus 
(1982), for S-degree latitude/longitude quadrangles centered at 37.5° N, 12.5° W 
(southern area) and 42.5° N, 12.5° W (northern area) have been plotted (Figures 5, 6, 
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Figure 3. ‘Time series of the six hourly values for the north/south component of the 
wind velocity (cm/s) from 1! February to 31 December 1986 at position 
39°N, 10°W: Day zero is | February 1986. 


and 7) and used to examine large scale features of the water properties in the study area. 
The upper water, from the surface down to about 125 m, includes the mixed layer and 
the seasonal thermocline (Figure 8). The laver between approximately 125 and 500 m, 
plotted in a temperature-salinity diagram (Figure 9), shows a linear relationship between 
temperature and salinitv with a slope of about 6.9°C per part per thousand (ppt). This 
is North Atlantic Central Water (Sverdrup e¢ al., 1942). The intermediate depths, from 
about 500 m and 2000 m, roughly the depth of the 35.1 ppt isohahine, are strongly in- 
fluenced by Mediterranean Water. The stronger influence of Mediterranean Water 1s 
located at 1000-1200 m and corresponds to a relative maximum of salinity (Figure 6), 
relatively high temperature (10.0-10.5°C, Figure 5), and potential density (o,) of about 
1027.75 kg/n? (Figure 7), as seen from climatology. Ambar and Howe (1979), have 
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Figure 4. Spectral densities of the north/south component of the wind velocity at 
39°N, 10°W from 1 February to 18 December 1986: Units are 


cnt s~? day for spectral density and day" for frequency. 


documented the separation of the Mediterranean outflow into two distinct cores, at 
about 7° W. Both cores flow westward at different depths, along the continental slope 
in the Gulf of Cadiz and turn northward around Cape Sao Vicente. At 38° N, off western 
Portugal, the upper core was observed with a relative temperature maximum of 
12.67°C at 660 m, and the lower core, recognized by a salinity maximum, was located 
further to the west at depths nearer 1300 m (Ambar and Howe, 1979). A third shallower 
vein of Mediterranean Water (distinct from the two main cores) was traced from the 


vicinity of the Strait of Gibraltar northward to about 38° 30’ N, off western Portugal, 


and followed the upper continental slope off the IP (Ambar, 1983). Below the 
Mediterranean Water, from 2000 m downward, North Atlantic Deep Water of relatively 
low stratification 1s present (Figure 7). Bottom waters with a temperature of about 
2.5°C and salinity of 34.9 ppt (Figures 5 and 6) are evolved forms of Antarctic Bottom 
Water (Tchernia, 1980). 

Comparing the climatological vertical profiles of im situ temperature, salinity, 
and potential density for the northern and southern areas, the following large scale fea- 
tures may be noted : 


¢ From the surface to about 2500 m, there is a substantial meridional gradient in the 
large scale temperature field (colder in the north, warmer in the south). From 2500 
m downward nearly uniform water temperatures correspond to North Atlantic 
Deep Water and Bottom Water. 


e Inthe upper 2500 m we note a substantial decrease of salinity from south to north, 
that 1s, a decrease of salinity with distance from the area of Mediterranean Water 
outflow (Figure 6). The meridional gradients of salinity and in situ temperature 
have the same vertical extent in the water column. Below 2500 m the salinity field 
looks uniform in the meridional sense. 


¢ The potential density field shows, in the upper 600 m, increasing values of o, from 
south to north, suggesting that the temperature effect overcomes the salinity effect. 
Bellow 600 m. one notes decreasing values of o, from south to north, especially at 
the levels of Mediterranean influence (1000-2000 m), indicating that salinity has an 
important control on density at these levels (Figure 7). 


3. The General Circulation 

Wooster er al. (1976), based on an analysis of mean monthly ship drift data, 
concluded that near the eastern boundary off Portugal there 1s maximum equatorward 
surface flow during the upwelling season. A seasonal reversal to predominantly 
northward winds during winter is accompanied by poleward alongshore surface flow. 
Current meter observations near the continental slope in 2450 m of water, 100 km west 
of Cape Sao Vicente, indicate northward flow at levels between 234 and 2046 m 
(Meincke et al. , 1975). In winter the countercurrent may reach the surface, but in other 
seasons, the poleward flow can have the configuration of a northward undercurrent be- 
low the upper 100 to 200 m, where the wind driven Portugal current flows southward 
(Fiuza, 1982). Moored current meter observations in 1100 m of water at 37° 24’ N, 9° 
25’W indicate a mean flow of 0.1 m/s towards the north over the period December 1982 
to May 1983 at 750 m; the meridional (v) component decrease with depth to almost 0 
at 200 m (Ambar, 1984; cited in Barton, 1989). Current meter data, in 1100 m of water, 
at 40° N, 9° 45’W, from July 1983 through January 1984, show a mean northward flow 
in excess of 0.3 m/s between 200 and 750 m (Ambar, 1985; cited in Barton, 1989). The 
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perature and depth have units of °C and meters, respectively. Data from 


Levitus (1982). 
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undercurrent off Portugal, deeper than 100-200 m, can flow northward and transport 


and water from a sub- 


water from the Mediterranean between 500 and 1500 m depth, 


tropical origin in the upper layers (Fiuza, 1982). 
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Figure 8. Climatological summer, winter, and annual temperature (°C) profiles in 
the upper 250 m of the water columm: for the 5-degree latitude/ long- 
tude quadrangle centered at 37.5° N, 12.5° W. Data from Levitus (1982). 
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Il. NUMERICAL MODEL 


A. MODEL DESCRIPTION 
I. Nfodel Equations 
The numerical model used in the experiments was developed by Haney (1985), 
modified by Batteen et al. (1989), and is a multi-level, primitive equation (PE) model. 


The hydrostatic, rigid lid, B-plane, and Boussinesq approximations are used in the model. 
The governing equations are as follows : 








du —1 op” 4 e7u 
7. ar + fv—A,Viut K,, Fr + 6,(u) {tj 
aol Ye (21 
dt Po oy a ” 62° 
z 2 “f 
ae eC US ects 
we | (Gee has \de [3] 
-H 
0 0 
p= | veae-+ | pedea: (4] 
Fis hanes 
= 
p = po(l — a(T — T)) i) 
aT oT 
Si 77 ANT + Ky 7 + Os + AT) . [6] 


In the above equations, t is time and (x,y,z) 1s a right-handed cartesian coordi- 
nate system with x pointing toward shore, y alongshore, and z upward. The corre- 
sponding velocity components are (u,v,w), T is temperature, p 1s density, and p’ 1s the 
departure of the pressure from the vertically averaged pressure. In equations (3) and (4), 
€ is a dummy variable of integration. Equation (4) includes the assumption that the 
depth-averaged pressure is a constant (assumed zero); that is, the barotropic mode 1s 


ignored. Equation (5) assumes that density is a function of temperature only. This as- 


Sumption may not be appropriate for the study region, because salinity also affects 
density over a substantial vertical extent due to the influence of the Mediterranean out- 
flow. The relatively important river runoff north of 41° N may also control the density 


near the coast. Future modeling studies of the region should incorporate salinity to as- 
Lameacl.S 
PoC Oz 








sess its effect on the circulation off the IP. In (6), Q,= is the heating (a 


function of depth) due to solar radiation, where 
S =S)(R e77! + (1 — R) e7/?2)_ [7] 


In formula (7), S, is the downward flux of solar radiation at the surface, R = 0.62 is the 
fraction of solar radiation absorbed in the upper few meters (z,=1.5 m), and 
(1 — R) = 0.38 1s the fraction that penetrates to somewhat deeper levels (z, = 20 m), as 
given by Paulson and Simpson (1977). The terms 6{u), 6,{v) and 6,{7) represent the 
vertical turbulent mixing of zonal and meridional momentum and heat, respectively, due 
to surface layer processes, by a dynamic adjustment mechanism. This adjustment main- 
tains the dynamic stability in the water column, in the sense of lighter water overlaying 
denser water. 


The boundary conditions at the top (z = 0) of the model ocean are : 


Ka = 0 [8] 
ON ea wees 
Km a = pp 9] 
) 
Ky ez —¥ Op [10] 
w=0 . [11] 


In (9), t is the alongshore component of the surface wind stress. In (10), Q, 1s the net 
upward flux of long wave radiation, sensible, and latent heat across the sea surface. 


The bottom (z = —//) boundary conditions are as follow : 


“A 


K, ae = Cy(u? + v*)!?(u cos y — v sin y) 2) 
ie = Gas =F yyy cosy +usin y) [13] 


one 


ae, 


= (0) [14] 


a. [15] 


In the above equations, Cp, is a bottom drag coefficient and y = 10° is the geostrophic 
inflow angle (Weatherly, 1972). Table 1 provides other symbols in the model equations, 


as well as values of constants used throughout this study. 


Table 1. VALUES OF CONSTANTS USED IN THE MODEL 


2 x da 
13x10 
Caensity ofair 


















JASE constant reference temperature 


: 2S erie density of air 


: 1.0276 gon cor density of sea water at 
2.0] <alieat As! thermal expansion coefficient 





iil 
Un 
No) 
— 
= 
= 
| 
i 


< AW ic 
K 
§. x10 on 
347 N 
aN 
Te 980. ems? ——SS~*dCnceeleration of gravy 


Ae 2. % 10M CHE biharmonic momentum diffusion coefh- 
cient 


2X AO Griese biharmonic heat diffusion coefficient 
Q.5 err? s7 vertical eddy viscosity 
Ke OUSFee is] vertical eddv conductivity 


2. Model Domain and Resolution 

The domain of the model is the rectangular region extending from approxi- 
mately 9° to 15° W and from 34.7° to 44.3° N (Figure 10), covering an area of 6 degrees 
of longitude by 9.6 degrees of latitude. The region extends approximately 500 km off- 
shore from the west coast of the IP, and it spans 1024 km in the alongshore direction, 
exceeding by about 1.3 degrees of latitude the northern limit of the coast (Cape 
Finisterre at 43° N), and by 2.3 degrees of latitude the southern limit (Cape Sado Vicente 
at 37° N). The horizontal grid spacing in the model is 8 km in the cross-shore direction 
and 16 km tn the alongshore direction, which should allow realistic, spatial resolution 
of mesoscale features in the Portugal Current System, since these features are expected 
to have typical wave lengths on the order of 100 km. Although there are variations in 
the coastline and ocean bathymetry (e.g., Cape Roca, and a few submarine canyons), 
these variations are omitted in the model to focus on the role of wind forcing in the 
generation of eddies. 

3. Finite Difference Scheme 

In the horizontal, a space-staggered B-scheme (Arakawa and Lamb, 1977; 
Batteen and Han, 1981) is used. There are 10 layers in the vertical, separated by con- 
stant z-levels placed at depths of 13, 46, 98, 182, 316, 529, 870, 1416, 2283, and 3656 m. 

4. Boundary Conditions 
The eastern boundary, representing the west coast of the IP, is modeled as a straight, 
vertical wall. A no-slip condition on the tangential velocity 1s invoked at the coastline. 
The northern, southern, and western boundaries are open using a modified version of the 
radiation boundary condition of Camerlengo and O’ Brien (1980). 

The 1624 km alongshore extent of the modeled coastline is not realistic, since 
the distance between Cape Sado Vicente and Cape Finisterre is approximately 667 km. 
However, the choice of this meridional extent allows wind forcing only along the corre- 
sponding section of the real coast. That is, the model uses both northern and southern 
bands of no wind forcing which, according to McCreary er al. (1987) and Batteen et al. 
(1989), will avoid unrealistic results for wind forced cases. The use of a band of wind 
forcing in the interior of the domain although somewhat artificial, allows for the propa- 
gation of coastal Kelvin waves, which can establish the alongshore pressure gradient 
field, resulting in a surface-trapped coastal jet and an approximately realistic undercur- 


rent (Batteen ef a/., 1989). 
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Figure 10. Mlodel domain area and ocean bathymetry: Depths in meters. 


5. Heat and Momentum Diffusion 
The model uses biharmonic lateral momentum and heat diffusion with the 


choice of coefficients listed in Table 1 . The choice of biharmonic, rather than 


Laplacian lateral diffusion, along with the appropriate coefficients, is justified since 
biharmonic diffusion is scale selective and acts predominantly on scales smaller than 
those of eddies (Holland, 1978), allowing mesoscale eddies to be generated via baroclinic 
and or barotropic instability processes (Holland and Batteen, 1986). 


B. INITIALIZATIONS AND SPECIFIC EXPERIMENTAL CONDITIONS 

The model can either be spun up from rest by a surface wind stress or heat flux, or 
it can be initialized with a specified current field. In this study the model is spun up from 
rest by a prescribed surface wind stress. 

1. The Initial Temperature Field 

Based on the mean climatological temperature profile (Figure 5) for the 5° 

latitude/longitude square centered at 37.5° N, 12.5° W, a simple approximation was used 
to initialize the temperature field in the model. This approximation consists of using 
exponentials for the upper 500 m and the lower 3300 m of the water column, with a 


linear variation between. The function fits the climatological profile and is given by : 


7(z)=9.84+ 7.7 ele (—500 m<z<0m) 
T(z) = 11.25 + 7.857 x 107°(z + 500),  (—1200 m < z < —500 m) [16] 


(2+1200) 
T(z)=2.5+8.2e' 600),  (—4500m<—1200 m). 


The resulting temperature profile (Figure 11) was introduced into the model to 
provide initial values at each of the 10 levels. The same initial temperature profile was 
used in each of the numerical experiments. The model computes density from temper- 
ature using equation (5). 

2. The Thermal Forcing 

The eastern boundary region off the west coast of the IP has, in the mean, a net 
heat gain from March through September (Bakun, 1987). This net gain occurs because 
of reduced latent heat flux and downward sensible heat flux in the coastal region, due 
to the presence of cold, upwelled water during summer. A surface thermal forcing 1s 
applied to simulate the real heat fluxes across the surface and avoid unrealistic low sur- 
face temperatures due to continuous upwelling forced by a constant equatorward, 
alongshore winds. 

The incoming solar radiation at the top of the atmosphere Q, , was specified to 
be 983.3 cal cnr? day~', as representative for June and July at latitude 39.5° N. The sum 


of the net long wave radiation, latent and sensible heat fluxes, Q, , was computed during 
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Figure 11. Initial temperature profile used in all experiments: Units for depth 


and temperature are meters and °C, respectively. 


the model experiments from standard bulk formulas (Haney et al., 1978), using April 
through September mean values of cloud cover, relative humidity, surface atmospheric 
pressure, the model simulated sea surface temperature field, the specified wind, and the 
air temperature. In the steady wind experiments, the air temperature was chosen initially 
to produce a zero net heat flux across the surface. Chmatological (April through Sep- 
teniber) values for cloud cover, relative humidity and surface atmospheric pressure of 
0.48, 86%, and 1018.3 mb, respectively, were taken from the climatological atlas by 
Esbenson and Kushnir (1981). Specified values for the wind speed and air temperature 


used for thermal budget computations in each of the experiments are indicated in 
Table 2. 
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Table 2. SPECIFIED VALUES OF WIND SPEED AND AIR TEMPERATURE 
FOR THE HEAT BUDGET OF THE EXPERIMENTS 


ee eee eee ee Eee 
Spatial mean wind speed 859 900 before day 20, 485 
(cm/s) 1425 after day 20 
Air temperature (°K) 289 289 before day 20, 290 |} 287 
after day 20 


3. Wind Forcing Experiments 











a. Experiment I (Large Negative Wind Stress Curl near the Coast) 
The wind field used in this experiment (Figure 12) consists of : 


e A strong northerly wind of 1426 cm/s (wind stress of about 3.5 dyne/cm’?) within 
48 km of the coast. 


e A sudden decrease of the alongshore wind stress to 754 cm/s (about 0.8 dyne/crm’) 
between 48 and 64 km off the coast. 


e Cyclonic wind stress curl, followed by anti-cyclonic curl, that approaches the 
climatological summer situation for the offshore part of the study region. 

In this experiment we add to the large scale climatological wind field, a 
strong negative wind stress curl near the coast which corresponds to the observations 
of McClain et al. (1986) along the west coast of Galicia. This feature, may be induced 
by orographic influences of the rias, and may not be representative of the entire west 
coast of the IP. However, we can consider this experiment to be process-oriented and 
focus on the resulting flow field, with the purpose of understanding the effect of these 
processes in eastern boundary current regions. 

Wind velocities are specified to be constant with latitude inside the forcing 
region. As previously discussed, there 1s no wind forcing in the northern and southern 
areas of the model domain. 

b. Experiment 2 (Poleward Wind Followed by Equatorward Wind) 

The purpose of this experiment is to simulate a “spring transition” and de- 
termine how the flow and temperature fields respond to such a transition. From day 0 
to day 20 the model is run with a steady, poleward wind of 900 cm/s (approximately 1.2 
dynes/cn?). From day 21 to the end of the model simulation (day 95), the model 1s forced 
with an equatorward steady wind of 1425 cm/s (approximately 3.5 dynes/cm’). Values 


of 1.2 dynes/cn?? and 3.5 dynes/cnv? are thought to be reasonable for peak values of strong 
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Figure 12. Wind forcing for experiment 1 


pulses of poleward and equatorward winds, respectively, for pre-upwelling and upwelling 
seasons. This wind forcing is used here to approximate, in a qualitative sense, the 
change from a situation of prevailing poleward winds to prevailing and intensified 
equatorward winds that may happen bv March. Note in Figure 3 and Table 4 (Chapter 
IV) that for the particular year of 1986, the prevailing poleward winds during February 
change to prevailing equatorward winds by the beginning of March, so that this exper- 
iment is not inconsistent with available observations. In other eastern boundary regions, 
such as off the coast of Oregon, “spring transitions” have been documented to corre- 
spond to significant changes in the therinal structure and velocity field over the coastal 


margin. 
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c. Experiment 3 (Forcing by a Time Series of the North[South Component of the 
Wind) 

This experiment consists of forcing the PE ocean model with the 
north/south component of the wind stress, computed from synoptic surface pressure 
analyses for the period from | February to 31 December 1986. This time series 
(Figure 3) was discussed in section I.B.1. This experiment constitutes the most realistic 
of the simulations, since it uses “realistic” winds derived from pressure observations in 


the study area. 


OR, 


Il. APPLICATION OF SPECTRAL ANALYSIS TO THE ESTIMATION 
OF THE DOMINANT SCALES AND ENERGY OF THE EDDY FIELD 


A. BACKGROUND 

In this section we describe the procedure and basis for a program which will be used 
to estimate the dominant wave lengths and kinetic energy of eddy fields that result from 
PE ocean model simulations. Spectral analysis is, in general, more suitable to estimate 
the spectral density function of random stationary processes. However, the method can 
also be used for deterministic processes, to pick out periodic components in the presence 
of noise (Chatfield, 1984). 

The PE model computes and stores, in memory, daily values for the fields of vari- 
ables (u,v, and T) at grid points identified by indexes (i,j,k). For variables u and v, both 
indexes 1 and j have the integer values 1,2,...,64. For variable T, the indexes 1 and } have 
integer values 1,2,...,65. Indexes i and j increase with eastward and northward directions, 
respectively. Indexes k= 1,2,...,10 (associated with the vertical) increase downward. In 
general (as Chapter IV will show), the resulting eddy fields tend to be generated and 
concentrated near the coast. The eddies tend to be anisotropic; perturbations develop 
preferentially along the y direction. For this reason, each spatial series is taken as a 
sample function corresponding to the values of u,v,I(j) (j= 1,2,...,.N), for a given (i,k); 
That is, discrete values along the northward direction (y) at level (k) and longitude index 
(i). Since in the model, the maximum j is 65 for temperature and 64 for the velocity 
components u and v, we used 64 points to compute the discrete Fourier transforms. 
Figure 13 shows, for purposes of illustration, a typical spatial series for the variable u 
along a north/south line near the eastern boundary , using data generated in experiment 
I. Notice in this particular meridional line that the zonal (u) velocity changes between 
positive and negative values with magnitudes and spatial scales that correspond to the 
expected intensity and wave lengths of mesoscale eddies. 

Before using the implemented program, one must consider several factors and decide 
which options to use. For instance, if there is relatively high variance at very low wave 
numbers and one wants to pick out a relative maximum in the spectrum at low wave 
numbers, one must consider the use of a prewhitening technique. This technique can 
help to pick out a relative maximum in the spectrum at higher wave numbers. I[t can 


also reduce the leakage of low wave number energy into higher wave number beams in 
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Figure 13. A representative spatial series of zonal (u) velocity (cm/s) along a 
north/south Hine in the model: The symbols, connected by straight 
lines, correspond to values of the variable (u) at grid points (i=S5, 


k= 1, j= 1,2,...,64) for experiment 1, on day 50. 


the presence of a very red spectrum (such as in our cases, which have a very red spec- 
trum for the high wave numbers). By differentiating the series, the spectrum can become 
more nearly white (flat) at the frequencies of interest, and smearing due to leakage of the 
energy in adjacent wave number beams can be reduced. A particular study case which 


tests the use of the different options is discussed later in this chapter. 


B. COMPUTATION OF THE SPECTRAL DENSITIES 

An individual spatial series of u,v, T(j) (j= 1,2,...,N), for fixed (i,k) will be referred to 
as g(i). The mean, B= Lol) is removed from each g(j) to get perturbations 
g'V)=e)—g. Ifa large trend shows up in the values (which usually does not happen 
in our experiments), g’(/) can be further detrended by subtracting the regression line 


values from each data point: 


aS 


Brew) = 26a) = (Og ray [17] 


The regression line coefficients, b, and m, are computed by the following formulas : 


N N 
> ve'() > 20) 
j=1 =| 


] 
m= (5) Baar narey ~ Wom a 


L\ oy 
b= W 2#0 ) = m(N = 1). [19] 


In the above formulas, dy is the grid spacing along the y direction, y = (j — l)éy, and 
N= 64 is the number of points in the sample to Fourier transform. 
To reduce the leakage problem due to the finite length of the series, one can use a 


cosine taper, u,, or a Hanning window, u,. The cosine taper window ts defined by: 


a 2} Y—1)5z so5-3NO 
u)=1~cos'| | i 





u(j=l, Ae<j<n- [20] 
w= 1-cos'| SSP? |, j2N-4. 
The Hanning window is defined by : 
u,(j) = | - cos!| x ont | La jaan [21] 


Then, g’new(/) = 2 od /UY), (= 1,2,....,N), where u(j) represents either window. 
The Fourier transform integral, applied to a finite record g’(v), O<y<JL,1s 


: . —l2nly 
xo=| ge dy, [22] 
0 


Where | 1s the y component of the wave number vector and ;=./—1. The above integral 


can be discretized as 
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. 
x) = Die’ exp] -1 EU NG - sy [sy (23) 


j=I 


/— | 
Now, / =1,2,...,.N is the index associated with wave numbers \ 7, ) and 





L =(N — !)dy is the length of the series. Wave number is defined as the inverse of the 
wave length. 


Dividing by dy, 


N 
= =) 2) exp| iar ¢—ne= |. fm 1,2,.0N, [24] 


J=1 


Summation on the right-hand side of (24) 1s computed by the discrete Fourier transform 


Fd 


oy 





subroutine for |= 1,2,....,N. Values are kept only for 1= 1,2,....N/2+1, that is, up 
to the Nyquist frequency. | 
The one sided spectral density for the series g(j), at discrete wave numbers ts given 


by 


KOK J 


C= ay? i= aaa [25] 


The subscript r denotes a particular sampling series along grid points j= 1,2,....,64 
and fixed indexes i,k. A superscript * denotes complex conjugation. 

At this point, to compensate for the loss of variance due to the application of a 
window, G,(/) is multiplied by a scale factor. The scale factors used are 2.0 and 1.1 for 
the Hanning window and cosine taper window, respectively, Which correspond to the 
ratio between the areas under the box car function (with a constant value equal to unity) 
along the record and under the window. 

In the case of stationary processes, to increase the confidence on the spectral esti- 
mates, the discretized Fourier transform integral is taken over a number of sampling 


series, NS, and the spectral estimates are averaged : 
NS 


C= DGD. [26] 


r=! 


Zi, 


Each sampling series contributes two degrees of freedom to the spectral estimates. Be- 
cause in the PE ocean model the variance of the field variables may change in the region 
Where the individual Fourier transforms are taken, the values G(l) represent how the 
spatially averaged variance 1s distributed in the wave number domain in the region of the 


model being considered. 


C. PREWHITENING 

The procedure for prewhitening consists of differentiating the series in the space 
domain and, after estimation of the one-sided spectral energy densities for the trans- 
formed series, using the transfer function to integrate in the wave number domain, which 
results in an estimate of the spectrum for the original data series. The application of this 
technique assumes that g(j) is a stationary process with a continuous power spectrum 
and that the transformation (differentiation) 1s stable, Ce) the variance of the output 
is finite. With these assumptions, the relation G(w) = wi holds (Chatfield, 1984). 


The differentiation is approximated by : 


7 £0) = 20") ee 
& ai {j = 1) im a see 2a 
by [27] 
2a ae 
After computing the spectrum G,,, the transfer function, w? = 2" Tp I 


1= 1,2,....N/2+1, is used to integrate the spectrum of the differentiated series to com- 


pensate for the prewhitening : 


a=. [28] 
w’ 


During the testing of the programs, several cases were run which involved changing 
the choices of the window (Hanning versus cosine taper) and experimenting with the use 
of prewhitening versus no prewhitening. Figures 14, 15, 16, and 17, show plots of 
spectral densities estimated for different combinations of prewhitening, no prewhitening, 
Hanning window, and cosine taper window. The sampling hnes g(j), }= 1,2,...,64 corre- 
spond in the model domain to indexes 1= 54,55,56,57, and k= 1,2. 

To help in the analysis of results regarding the conservation of variance for the dif- 
ferent cases, Table 3 shows the values for the ratios of variances associated with u, v, 
and T, as computed from integrating the averaged spectrum, G(l), and the averaged 


Variances as computed directly from the individual samples, g’(/). 
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Figure 14. Spectral densities with prewhitening and Hanning window, for zonal ve- 
locity and temperature for experiment [ on day 50: Units are 


(cin/s)? km and (°C)?km for the spectral densities of velocity and tem- 


perature, respectively. 
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Spectral densities without prewhitening but with Hanning window, for 


Figure 15. 
Units are 


zonal velocity and temperature for experiment I on day 50: 
(crt/s)? km and (°C)? Art for the spectral densities of velocity and tem- 


perature, respectively. 
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Spectral densities without prewhitening but with cosine taper window, for 


Figure 16. 
Units are 


zonal velocity and temperature for experiment 1 on day 50: 
(cm/s)? km and (°C)? km for the spectral densities of velocity and tem- 


perature, respectively. 
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Day 50: Region (= 54 ,57; J=1,64; K=1,2) 
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Figure 17. Spectral densities with prewhitening and cosine taper window, for zonal 
velocity and temperature for experiment 1 on day 50: Units are 
(cri/s)? Am and (°C)? Am for the spectral densities of velocity and tem- 


perature, respectively. 
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Table 3. RATIO OF VARIANCES AS COMPUTED FROM THE_ INTE- 
GRATION OF THE AVERAGED SPECTRUM AND AS COMPUTED 
DIRECTLY FROM THE SPATIAL SERIES 


 -_ ar a ee 
Prewhitening and Hanning 
No Prewhitening and Hanning eae 
No Prewhitening and Cosine Taper 0 i o.09 
Re 










V 


Prewhitening and Cosine Taper 





Looking at the summary of results, the cosine taper window does a better job (compared 
with the Hanning window) in conserving variance in both the prewhitening and no 
prewhitening cases. Notice that the average value of the ratios shown in Table 3 1s 
closer to unity for the cases using the cosine taper window. In general, the ratios shown 
in Table 3 are relatively close to unity, except for the variable u in the cases associated 
with the first two rows of Table 3 (in which a Hanning window was used). In those 
cases, the values for the ratios are 1.36 and 1.41, and are probably due to overcom- 
pensation of variance by applying the scale factor of 2.0 for a Hanning window. From 
these results it follows that a cosine taper window and no prewhitening should give 
better estimates for the mean eddy kinetic energy computed by integrating the spectral 
densities. 

From the plots of spectral densities for variables u and T (Figures 14, 15, 16, and 


17) one notes the following : 


¢ Comparing cases without prewhitening and with prewhitening using the cosine ta- 
per (Figures 16 and 17), it is seen that the shape of the spectrum is about the same; 
prewhitening does not seem to affect the results in a considerable manner. 


¢ Comparing the cases of no prewhitening using the Hanning window and the cosine 
taper window (Figures 15 and 16), notice that the Hanning window yields a 
smoother spectrum at the higher frequencies due to larger suppression of the side 
lobes in the Fourier transform. 


¢ Comparing the cases of the Hanning window with prewhitening and without 
prewhitening (Figures 14 and 15), there is not much difference in the shape of the 
spectrum, except for very low wave numbers. For the very low wave numbers, the 
prewhitening appears to smooth the spectrum for the variable u but not for T. 
lTowever, the cosine taper does a better job in separating two peaks at very close 
frequencies (resulting in better resolution) due to a narrower mam lobe of the 
Fourier transform of the cosine taper window compared with that of the Hanning 
window. 
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From the above considerations, it is suggested that if one desires to pick out the relative 
maximum from the spectrum, one may get better results using a Hanning window either 
with prewhitening or without prewhitening, since it does a better job in reducing leakage, 
unless the relative maxima are very close in wave number (the use of a logarithmic scale 
helps to show more detail in the spectrum over a wider range of wave numbers, com- 


pared to a linear scale). 


D. EDDY KINETIC ENERGY 

Having computed the spectral densities for fluctuations u’, v’, and 7’, for a given 
region in the model, and for a given day, we may look towards the computation of the 
daily sequence of the spatial mean eddy kinetic energy. To do this, take 
= © G.,(081 to approximate KE’,,. The value, 6/= —, is the wave number resol- 
ution, and KE’, represents the contributions of u’,v’ to the mean eddy kinetic energy. 
Then KE’,,,,. = 02 + o? represents the spatial mean (in the region defined by 1,),k) eddv 
kinetic energy associated with both wu’ and v’ (an example is shown in Figure 18). If one 
desires the eddy kinetic energy for interpretation of experimental results, it 1s more ac- 
curate to compute the energy directly from the fields of u and v, rather than by inte- 
grating the spectral densities. 

A daily sequence of the spectral estimates for u’ is shown in Figure 19. The choice 
of a linear rather than a logarithmic scale allows easier visualization of the dominant 
peaks. The combination of Figure 18 and Figure 19 shows very clearly how the spatially 


averaged eddv kinetic energy and dominant wave lengths can change with time for a 


given region of the model, for a particular simulation. 
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EDDY KINETIC ENERGY 
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Time sequence of eddy kinetic energy computed by adding the integrated, 


spatially averaged spectral densities of zonal (u) and meridional (v) com- 
Units are cn? s-? for eddy kinetic energy and 


Figure 18. 


ponents of the velocity: 


days for time. 
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Figure 19. Time sequence of the spectral densities of variable (u), constructed by 
averaging the spectral densities of alongshore lines of grid points for ex- 
periment 1, from day 30 to 100: Units are (cm/s)? kim for the spectral 


densities of the (u) velocity and days for time. 
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IV. RESULTS OF WIND FORCING EXPERIMENTS 


A. STEADY WIND FORCING EXPERIMENT 
1. Experiment 1 (Large Negative Wind Stress Curl near the Coast) 

A similar experiment was conducted by McClain er al. (1986) to study the flow 
field forced by of a negative wind stress curl observed off the west coast of Galicia. The 
diagrams of the applied wind forcing and the resulting current velocity field and 
interfacial displacements for the McClain et al. (1986) experiment, are shown in Figures 
20 and 21. The southward flow adjacent to the coast was driven by the southward 
coastal wind and the northward flow offshore was associated with the negative curl of 
the wind stress (McClain et al., 1986). This same type of structure in the wind field was 
used in experiment | for the study region off the west coast of the IP. 

In experiment | the wind forcing (Figure 12) consisted of : a large negative wind 
Stress within 48 km of the coast; a large negative wind stress curl in the band between 
48 and 64 km from the coast; a positive wind stress cur! between 64 and 184 km; and 
offshore of 184 km, a negative curl, the latter structure based on the climatological wind 
Stress curl for the IP during the upwelling season (Bakun, 1987). These winds change 
only in the cross-shore direction and are steady in time. As in all other experiments, only 
the alongshore component of the wind 1s used. 

a. Development of the Coastal Jet and Offshore Countercurrent 

On day 5, due to the equatorward wind stress, a surface westward flow 
(Figure 22 (left)) results, due to Ekman transport away from the coast, as expected. The 
largest negative values (magnitudes over 10 cm/s) for the flow are seen in a meridional! 
band with its axis located approximately 60 km from the coast. Also seen in Figure 22 
(right) is the surface dynanuic height field (relative to 2400 m), which shows a ridge, with 
its axis extending alongshore approximately 55 km from the coast. A cross-section taken 
in the middle of the domain (Figure 23), shows that the upper ocean isotherms slope 
downward toward the coast at distances between approximately 82 and 57 km from the 
coast and slope upward within 57 km of the coast. The opposite slopes for the isotherms 

Ou 


fall in regions where the gradient ( oe ) at the surface has opposite signs (compare Fig- 


ure 23 and Figure 22 (left)). The change in sign of the zonal gradients of the surface 





dynamic height field, of the zonal velocity, and of the upper ocean temperature field, 


corresponds to the region where the wind stress curl 1s negative (1.e., between 48 km and 
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Figure 20. Numerical model basin geometry and the definition of the southward 


wind field applied in the simulation of McClain et al. (1986). 


64 km from the coast). The upwelled isotherms near the coast are consistent with a di- 
vergence of the flow in the Ekman layer under the influence of the equatorward wind 
with a compensating vertical motion to replace the surface water. Farther offshore, with 
the decreasing Velocities of the equatorward wind, the offshore Ekman transport is re- 
duced resulting in surface layer convergence and downwelling consistent with the pattern 
of isotherms seen in Figure 23. Because the surface dynamic height field is computed 
by vertically integrating the specific volume anomaly (which is directly proportional to 
the temperature), the ridge shown in Figure 22 (right) coincides with the location where 
upper ocean isotherms have a maximum downward displacement. Further, the region 
of imposed negative wind stress curl (distances from the coast between 48 and 64 kin) 
coincides with the locations of both the ridge of the surface pressure and the change in 
sign of the slopes of the upper ocean isotherms. The pressure gradient force, due to the 
surface slope and horizontal gradient of temperature/density, changes in sign across this 


region of negative wind stress curl so that there results an equatorward surface jet within 
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Figure 21. Thie flow field and pycnocline displacement resulting from the wind forc- 
ing shown in Figure 20: Dashed lines are lines of constant displace- 
ment of the interface, and the shaded area defines the area of strong 
wind stress curl. The maximum wind stress is 3 dynes/cm (from 
McClain et al., 1986). 


approximately 55-60 km from the coast and a poleward surface jet offshore. This ts 
clearly identified 1n a zonal vertical section of the alongshore averaged v component of 
the velocity (Figure 24), which shows the separation between jets at about 57 km from 
the coast. The equatorward jet has an alongshore averaged maximum of about 30 cm/s 
at 10 km from the coast, and a depth of about 320 m. The poleward jet has a maximum 
alongshore averaged velocity of about 15-20 cm/s at 67 km from the coast and a vertical 
extent of approximately 220 m. Results of this experiment on day 5 are qualitatively 
consistent with the results of McClain et al. (1986), in which the currents forced under 
a similar wind field reached a maximum southward speed of 15 cm/s shoreward of the 


wind shear zone and a maximum northward speed of 14 cm/s offshore of the zone. 
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Figure 22. Surface isopleths of (left) zonal (u) velocity (cm/s) and (right) dynamic 
height (cm) relative to 2400 m for experiment | at day 5: The contour 
interval 1s 5.0 cm/s for velocity and 2.0 cm for dynamic height. Dashed 


lines denote offshore velocities and negative values relative to 2400 m. 


By day 20, the coastal equatorward surface jet, within 60 km of the coast, 
has an intensified maximum alongshore averaged speed of about 60 cm/s and has 
deepened its vertical extent to about 860 m depth. The poleward jet also has an in- 
creased maximum alongshore averaged speed of about 35-40 cm/s, and has increased Its 


cross-sectional area by deepening to 560 m depth, and widening to about 123 km from 
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Figure 23. Vertical cross-shore section of temperature (°C) at y=512 km for ex- 


periment I at day 5: Contour interval is 1.0°C. 


the coast. A new, although weak (approximately 5 cm/s), feature in the system is a 
poleward undercurrent, which shows up adjacent to the coast between approximately 
200 and 500 m (Figure 25). The equatorward surface flow and poleward undercurrent, 
which result from forcing by an equatorward wind stress in the coastal region, are con- 
sistent with the results of McCreary et a/. (1987) and Batteen et al. (1989). 
b. The Generation of Eddies 

Superimposed on the dominant poleward and equatorward surface jets 
(Figure 26 (left)) are eddies, as seen in the surface zonal current field (Figure 26 (right)) 
at approximately y= 735 and y= 191 km. These perturbations are confined to within 120 
km of the coast in the same region as the poleward and equatorward jets as shown in 
Figure 26(left). 
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Figure 24. Vertical cross-shore section of alongshore-averaged meridional (v) veloc- 
ity (cm/s) for experiment 1 at day 5: Contour interval is 5.0 cm/s. 


Dashed lines denote equatorward flow. 


Following Batteen er al. (1989), we investigate the necessary conditions for 
instability by examining the distribution of potential vorticity, q, 11 ume-averaged 
cross-sections of the coastal flow. Following Watts (1983), the potential vorticity, q, 


may be approximated by the expression: 











OT oT ov 
q= (ft ¢) G2 ook ax Az ’ [29] 
where pat Se A cross-section of the time-averaged (days 21-30) potential 


vorticity at y= 768 km, is shown in Figure 27. The choice of davs 21 through 30 for the 


time averaging corresponds to the period in which the instabilities developed. Away 
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Figure 25. Vertical cross-shore section of alongshore-averaged meridional (v) veloc- 
ity (cm/s) for experiment 1 at day 20: Contour interval 1s 5.0 cm/s. 


Dashed contours denote equatorward velocities. 


from the coastal region, there is a tendency for potential vorticity to be uniform along 
isothermal surfaces and to change vertically, consistent with the temperature 
Stratification (not shown). A relative maximum of potential vorticity 1s present at a 
depth of about 85 m, corresponding to the seasonal thermocline (1.e., the maximum 
vertical gradient of temperature). In the upper ocean, within 128 km of the coast, the 
distribution of isolines of potential vorticity is complicated by large horizontal gradients 
of temperature and velocity shears. The cross stream derivative of potential vorticity 
(calculated by computing the horizontal derivative (s*) and then multiplying by one 
grid length (dx) ) changes sign in the upper ocean (Figure 28), at a distance of 54 km 


from the coast, corresponding well to the axis where the anticyclonic eddies are gener- 
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Figure 26. Surface isopleths of (left) meridional (v) velocity (cm/s) and (riglit) zonal 
(u) velocity (cm/s) for experiment 1 at day 20: Contour interval 1s 10.0 


cm/s. Dashed lines denote equatorward and offshore velocities. 


ated in the experiment. Because there is a sign change in the zonal direction at this 
distance froin the coast, which satisfies a necessary condition for barotropic instability, 
barotropic instability is likely the main mechanism for the generation of the anticyclones 
in experiment 1. Note however, that there is also a sign change in the vertical direction, 


which satisfies a necessary condition for baroclinic instability. As a result, mixed 
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(barotropic and baroclinic) instability could lead to the generation of eddies in this ex- 
periment. 
c. The Eddy Fields 

By day 40, the surface dynamic height field (Figure 29 (left)) shows the ex- 
istence of four, well developed anticyclonic eddies with center positions at y ~ 252, 438, 
604, and 755 km. Because the centers of the eddies nearly coincide with the zone of large 
negative wind stress curl, the four anticyclonic circulations are likely generated due to 
the anticyclonic shear of the wind in this region. The three northernmost eddies have 
circular patterns, with diameters of approximately 147 km, while the southern eddy has 
an elliptical pattern, with its major axis oriented northeast/southwest. All of the eddies 
have warm cores (consistent with warm core rings) as seen in Figure 29 (right) and in 
Figure 30 for the eddy located at y ~ 624 km. The near-surface isotherms associated 
with the eddies show considerable downward displacements. For instance, as shown in 
Figure 30, the 17°C isotherm surrounds the eddy at the surface, and is displaced down- 
ward to 186m. The surface temperature pattern (Figure 29 (right)) has four cold water 
filaments corresponding to places where the circulation of each of the four anticyclones 
advects the upwelled coastal water offshore. The southern filament is the largest, 
advecting relatively cold water (below 15.0°C) as far as 185 km offshore. 

The spectral densities of the fluctuations of u,v, and T associated with the 
alongshore component wave number were computed by sampling the fields at the model 
grid points (upper level, O< y < 1008 km, distances from the coast between 32 km and 
256 kin) corresponding to the region with the most eddy activity and are shown in 
Figure 31. Notice that the peaks of the spectra for u,v, and T show the same alongshore 
spatial scales. The dominant meridional wave number is approximately 0.0065 km-', 
which corresponds to a wave length of 154 km. This agrees with the scaies of the eddies 
seen in the surface dynamic height field (Figure 29(left)). 

By day 60, the two northern and the two southern eddies have 
merged/coalesced into one northern and one southern eddy (Figure 32 (left)). The 
dominant alongshore wave number of the eddies (Figure 33) 1s 0.006 Am! (alongshore 
wave length of 166.7 km). The resultant eddies continue to enlarge, as shown for day 
90 (Figure 32 (right)). This expansion is also seen in the time series of daily spectral 
densities of the u component of the velocity field (Figure 34). 

Note that the southern eddy extends farther offshore than the northern 
eddy. This can be explained by the f-effect. The two anticyclones propagate westward 


at estimated rates of about 1.4 km/day and 2.3 km/dav for the northern and southern 
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Figure 27. Vertical cross-section of potential vorticity (°C /(m s)) for the time- 
averaged days 21-30 of experiment 1: Values are scaled bv 10° . Con- 
tour interval is 0.2°C / (ms). The vertical cross-section was taken at 
y= 768 km. 


anticyclones, respectively. The faster rate of the southern eddy is consistent with a 
B-induced motion, 1.e., 8 is larger at lower latitudes. 

By day 90, the size of the eddies has stabilized. The diameter of the eddies, 
defined by the 10 cm dynamic height contour, is approximately 166 km (Figure 32 
(right)), while the dominant meridional scale of the eddy field is around 300 km, which 
corresponds to an alongshore wave number of 0.0033 kni', as shown in the spectral 
density distribution (Figure 35). A vertical zonal section in the middle of the domain 
shows the thermal structure of the northern warm core ring as a relatively large down- 


ward displacement of the upper ocean isotherms (Figure 36). For instance over a zonal 
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Figure 28. ‘Vertical cross-section of the cross-stream derivative of potential vorticity 
multiplied by the zonal grid size (°C / (im s)) for the time-averaged days 
21-30 of experiment 1: Values are scaled by 10°. Contour interval is 
0.1°C / (ms). Dashed contours denote negative values. The vertical 


cross-section was taken at y= 768 km. 


distance of 74 km, the 17°C isotherm bends down from the surface to 223 m and then 
bends up to the surface. This thermal structure represents a considerable amount of 
available potential energy for converting to kinetic energy associated with this particular 
eddy. [The vertical structure of the v component of the velocity at the same alongshore 
location (Figure 37) is consistent with the thermal structure (via the thermal wind re- 
lation), and shows a poleward flow west of the center position of the eddy (located at 
150 km from the coast) and an equatorward flow on the eastern side. In addition, there 


is a coastal surface jet coincident with the upward slope of isotherms in the nearshore 
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Figure 29. Surface isopleths of (left) dynamic height (cm) relative to 2400 m and 
(right) temperature (°C) for experiment 1 at day 40: Contour interval! 
is 5.0 cm for dynamic height and 0.5°C for temperature. Dashed con- 


tours denote negative values relative to 2400 m. 


region. This jet, which has a maximum surface velocity of about 60 cm’s, interacts with 
the offshore southward flow of the anticyclonic eddy circulation (at approximately 67 
km offshore) resulting in a deep southward flow at this position (v 1s larger than 10 cm/s 
as deep as 1300 m). By dav 90, a plume of cold water (temperatures less than 16°C ex- 


tending approximately 250 km offshore) lies between the two anticyclones at y= 416 km 
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Figure 30. Vertical cross-shore section of temperature (°C) at y = 624 kin for ex- 


periment | at day 40: Contour interval is 1.0°C. 


(Figure 38), and appears to have been advected by the southward branch of the north- 
ern eddy. A zonal section of the temperature field across the cold filament shows 
sinusoidal-like isotherms superimposed on a general upward sloping towards the coast 
(Figure 39). A cold dome (with an upward displacement of the 15°C isotherm at about 
202 km offshore) is associated with a cyclonic meander that is evident in the surface 
dynamic height field (Figure 32 (right)). A comparison between the spectral densities 
of the u component fluctuations for days 60 and 90 (Figures 33 and 35) shows that the 
dominant alongshore component wave number has changed with time towards lower 
values, resulting in an increase in the meridional scale of the eddies. 

From day 90 to the end of the experiment (day 154), the two anticyclones 


continue to propagate westward at rates of approximately 2.2 and 3.1 km/day, respec- 


49 


Doy 40: Region (= 33 ,64 J=1, 64; K =i) 


10000000 


> 
ze 
Uv} 
c 
D 
a) 
re) 
Ps 
5) 
@ 
Q 
Vv? 





Figure 31. 
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Spectral densities for the alongshore wave number of zonal (u) velocity, 
meridional (v) velocity, and temperature (I) for experiment 1 at day 
40: The spectral density estimates were obtained from the spatial se- 
ries of variables u,v, and T at grid points along the meridional lines 
defined by indices (1= 33,34,...,61; j= 1,2,..,64; k= 1). Units of spectral 


densities are cnr s-? kt for u and v, and (°C)? kon for T. 
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Figure 32. Surface isopleths of dynamic height (cm) relative to 2400 m for day 60 
(left) and day 90 (right) for experiment 1: Contour interval 1s 5.0 cm. 


Dashed lines denote negative values. 


tively, for the northern and southern eddies. This westward propagation of the eddies 
can be seen by comparing the surface dvnamic height field in Figure 32 (right) with 
Figure 40. By day 150, the two large eddies have propagated farther offshore and two 
new anticyclones appear in the same region of negative wind stress curl as before, at 
approximately y= 768 and y= 384 km. The thermal structure (not shown) of these new 


warm core rings is consistent with the velocity field (via the thermal wind relation). 
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Figure 33. Spectral densities for the alongshore wave number of zonal (u) velocity 
and temperature (T) for experiment 1! at day 60: The spectral density 
estimates Were obtained from the spatial series of variables u and T at 
grid points along the meridional lines defined by indices (1= 33,...,61; 
j= 1,2,...,64; k= 1,2,3). Units of spectral densities are crm? s-? kit for u 


and (°C)? km for T. 
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Figure 34. Spectral density versus time of the alongshore wave number of zonal (u) 


velocity for experiment I from day 60 to day 95: The spectral density 
estimates were obtained from the spatial series of (u) at grid points 
along meridional lines defined by indices (1=6,7,...,64; j= 1,2,...,64; 


k=1). Units are cn? s-? kom. 


B. TRANSIENT/TIME DEPENDENT WIND FORCING EXPERIMENTS 


1. Experiment 2 (Poleward Wind Followed by Equatorward Wind) 


a. Spring Transition 
Forcing in this experiment consisted of a spatially uniform poleward wind 
stress of 1.2 dynes/cm? during the first 20 days followed by spatially uniform equatorward 
wind stress of 3.5 dynes/cn? to the end of the experiment (day 105). This wind forcing 
is used to qualitatively represent the transition between a pre-upwelling season and the 


upWelling season that may happen in March or April. The values of wind stress, rather 
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Figure 35. Spectral densities for the alongshore wave number of zonal (u) velocity 
and temperature (T) for experiment | at day 90: The spectral density 
estimates were obtained from the spatial series of variables of u and T 
at grid points along the meridional lines defined by indices (1= 33,...,61; 


j= 1,2,...,64; k=1,2,3). Units of spectral densities are ci? s-? Am for u 


and ( sGo2screl Ohms 


S4 


T PEED at Day SO°at y—1S12 km 


West—east distance (km) 
320 256 192 


=> 
£ 
= 
es 
v 
(an) 


contour interval = 1.0 





Figure 36. Vertical cross-shore section of temperature (°C) at y=512 km for ex- 


periment 1 at day 90: Contour interval ts 1.0°C. 


larger than the mean climatological values for the region, correspond to peak values of 
strong wind events in the region for pre-upwelling and upwelling seasons. 

Under the forcing of the poleward winds, as expected, the surface flow re- 
sponds with onshore Ekman transport (Figure 41 (left)), convergence in the surface 
laver and downwelling, as seen by the tilting of the 17, 16, and 15°C isotherms in the 
upper 200 m, within approximately 20 km of the coast (Figure 42). On dav 20 the sur- 
face dynamic heights relative to 2400 m show increasing values towards the coast (not 
shown) and a coastal surface poleward jet consistently shows up within about 40 km 
from the coast, with increasing velocities from south to north ( Figure 41(nght)). The 
alongshore averaged meridional (v) velocities ( Figure 43) show that the vertical extent 


of the surface poleward jet 1s about 150 m, that the core is at the surface at about 8 km 
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Figure 37. Vertical cross-shore section of meridional (v) velocity (cm/s) at y = 512 


km for experiment | at day 90: Contour interval is 10.0 cm/s. 


offshore, and that there is a Weak equatorward undercurrent with an alongshore aver- 
aged velocity of around 2 cm/s at depths between roughly 200 and 520 m, corresponding 
to the portion of the section where the 14, 13, and 12°C isotherms slope slightly upward 
towards the coast (Figure 42). While the equatorward undercurrent is not commonly 
observed near the coasts of eastern boundary regions, with the exception of the Leeuwin 
current (Thompson, 1984) the surface poleward current under poleward winds is con- 
sistent with observations on the shelves. For instance, according to Hluyer (1983), over 
the shelf off Oregon during the winter (prevailing poleward winds), there is no mean 
southward surface flow and no mean vertical shear. The simulated poleward current in 
the upper 150 m partially corresponds to those observations. However, the vertical shear 


of the poleward current in the model does not correspond as well to Hluvyer’s (1983) ob- 
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Figure 38. Surface isotherms (°C) for experiment 1 at day 90: Contour interval 
sQaer C. 


servations. This difference may be caused by the exaggerated downwelling of the 
isotherms in the model for the coastal region, leading to large horizontal density gradi- 
ents in the offshore direction and associated vertical shear in the meridional (v) veloci- 
ties. The core of the poleward jet is placed approximately 8 km offshore in the position 
Where the 17°C isotherm slopes downward from the surface. On actual continental 


Shelves, the location of the core of the poleward jet may be modified by seasonal river 
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Figure 39. Vertical cross-shore section of temperature (°C) at y=416 km for ex- 


periment ! at day 90: Contour interval is 1.0°C. 


discharges that introduce changes in the honzontal stratification of the near coastal 
waters, but such effects are not addressed in this study. 

After the transition to equatorward winds on day 22, a vertical zonal section 
of temperature in the middle of the domain (Figure 44) shows that the previous down- 
welled isotherms now tilt upward, particularly near the coast. However, the 15 and 
16°C isotherms, at distances from the coast between 26 and 10 km, still show the 
downward tiling which remains from the previous forcing by poleward winds. By day 
22, the poleward surface jet evident on day 20 is not apparent and has been replaced by 
a relatively deep (vertical extent of about 350 m at the coast), equatorward coastal cur- 
rent (Figure 45). The honzontal maps of zonal (u) velocity at depths of 20 and 70 in 


(Figure 46) show that the surface Ekman laver may get as deep as about 70 1m in the 
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Figure 40. Surface contours of dynamic height (cm) relative to 2490 m for exper- 
iment 1 at day 150: Contour interval is 5.0 cm. Dashed contours de- 


note negative values. 


model. At 20 m, relatively large values of offshore zonal (u) velocity (Figure 46 (left)) 
are consistent with a typical Ekman laver response to equatorward winds. Below the 
Ekman layer (depths of typically 20-30 m), observations on the shelves in eastern 
boundary regions have shown onshore transport in the so called “geostrophic interior 
regime’. [he model seems to reproduce a relatively deep Ekman layer, because as dem- 


onstrated in Figure 46 (right), the transition between the two regimes may happen at 
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Figure 4]. Surface isopleths of (left) zonal (u) velocity (cim/s) and (right) meridional 
(v) velocity (cm/s) for experiment 2 at day 20: Contour interval is 2.0 
cm’s for (u) and 5.0 cm/s for (v). Dashed lines denote offshore and 


equatorward velocities. 


about 70 m depth. Comparing the results described for days 20 and 22 of this exper- 
iment, one may conclude that the model coastal isotherms and velocities respond rela- 
tively quickly (one or two days) to a sudden change from moderate poleward winds to 


strong equatorward winds. 
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Figure 42. _— Vertical cross-shore section of temperature (°C) at y=512 km for ex- 


periment 2 at day 20: Contour interval is 1.0°C. 


On day 30, after 10 days of equatorward winds, the 16 and 15°C upwelled 
isotherms have surfaced at distances of about 30 and 13.5 km from the coast, respec- 
tively (Figure 47). The surface equatorward jet, with its core located approximately in 
the position where the 15°C isotherm has surfaced, shows a large vertical shear in the 
alongshore (v) velocity (Figure 48) that is associated with the frontal-like structure in 
the upper 70 m of the nearshore water column. In this region, temperature increases in 
the offshore direction by 2.0°C in 30 km (Figure 47). A poleward undercurrent, with 
an averaged meridional (v) velocity greater than 5 cm/s, also exists nearshore between 
320 and 450 m (Figure 48), where the 12°C isotherm bends down towards the coast 
(Figure 47). By day 30, the system consists of an equatorward surface jet, located within 


40 km of the coast, and a poleward undercurrent at depths greater than 186 m, right at 
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Figure 43. Vertical cross-sliore section of alongshore-averaged meridional (v) ve- 


locity (cm/s) for experiment 2 at day 20: Contour interval is 2.0 cm/s. 


the coast. At this point we may establish a parallel with observations over the shelf off 
Oregon that document a “spring transition” under a single upwelling event that began 
on 25 March and continued intermuttently through 31 March, 1975. According to Huver 
(1983), this particular upwelling event “caused the transition from the typical winter re- 
gime with level isopycnals and no mean vertical shear to the typical spring/sununer regime 
with sloping isopycnals, a mean southward surface current and a Strong niean vertical 
shear”. We note from our experiment the following differences relative to the tvpical 
winter situation off Oregon : 


e Before the turning to equatornward winds, the isotherms/isopycnals were sloping 
downward near the coast as opposed to being level. 
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Figure 44. Vertical cross-shore section of temperature (°C) at y=512 km for ex- 


periment 2 at day 22: Contour interval is 1.0°C. 


e Before the turning to equatorward winds, there was a vertical shear in the poleward 
current nearshore, associated with downwelling, as opposed to no mean vertical 
shear. 

Nevertheless, starting from slightly different conditions, after 10 days of 14.2 mvs, 
equatorward winds (in the observations off Oregon, the southward upwelling favorable 
Winds were intermittent with a peak velocity of about 12.7 m/s), the model thermal and 
velocity fields show a change that is consistent with the spring transition observations 
off Oregon; by day 30, there are upward sloping isotherms/isopycnals, a mean southward 
surface current, and a strong mean vertical shear. For the study region off the coast of 


Portugal, there is a lack of observational studies to document the seasonal transition 
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Figure 45. Vertical cross-shore section of alongshore-averaged meridional (¥) ve- 
locity (cm/s) for experiment 2 at day 22: Contour interval is 5.0 cms. 


Dashed lines denote equatorward velocities. 


that may occur in March or April and thus we are not sure if a single strong wind event 
will mark the beginning of the upwelling season. 

With time, under the forcing of the steady equatorward winds (compare 
Figures 47 and 48 with Figures 49 and SO, for days 30 and 50) we note the following 
changes: 


e The surfaced isotherms are displaced farther offshore and the core of the 
equatorward jet 1s also displaced offshore, apparently tracking the 14.5°C isotherm. 


e The jet widens as far as 64 km from the coast and its maximum alongshore 
meridional (v) velocity weakens from ~ 35 cm/s on day 30 to ~ 28 cm,s on day 50. 
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Figure 46. —_Isopleths of zonal (u) velocity (cm/s) at 20 m depth (left) and 70 m depth 
(right) for experiment 2 at day 22: Contour interval is 2.0 cm/s. 


Dashed lines denote offshore velocities. 


e The depth where the southward current changes to a poleward current 1s seen to 
shoal somewhat within approximately 20 km of the coast. For instance, on day 30, 
poleward flow existed at depths greater than about 185 m, while on day 50, the 
poleward flow was below about 120 m. 


e A relatively weak surface poleward countercurrent exists offshore of the 
equatorward surface jet on day SO. 
At this point, the results of the model suggest that as one progresses through the up- 


welling season, the surface layer front and the core of the equatorward surface jet may 
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Figure 47. =‘ Vertical cross-shore section of temperature (°C) at y= 512 km for ex- 


periment 2 at day 30: Contour interval is 1.0°C. 


tvpically be displaced offshore (as compared to earlier in the season) and that a com- 
pensating poleward undercurrent may reach depths that would correspond to the lower 
portion of a ccntinental shelf. According to Huver (1983), “the dominant features are a 
southward coastal jet at the surface and a poleward undercurrent along the bottom {and} 
the southward coastal jet has its maximum (about 25 cm/s) at 15-20 km from shore [and] 
the poleward undercurrent has a velocity of about 5 cm/s’. Current meter observations 
taken off the coast of Portugal at latitude 38°06’N, from 12 June to 7 July of 1980, 
during the peak of the upwelling season, had mean meridional (v) velocities of -1.69 cm/s 
at 15 m depth, 4.6 km offshore (bottoin depth of 50 m) and -2.08 cm’s at a depth of 17 
m, 21 km offshore (near the shelf break and bottom depth of 150 m; Fiuza, 1984).Thus 


the model results on day SO (after 30 davs of equatorward winds) seem consistent ob- 
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Figure 48. Vertical cross-shore section of alongshore-averaged meridional (¥) ve- 
locity (cm/s) for experiment 2 at day 30: Contour interval! is 5.0 cm/s. 


Dashed lines denote equatorward velocities. 


servations off Oregon, by showing the equatorward surface jet with the core at a distance 
of about 27 km, a maximum velocity of about 28 cm,s, strong vertical shear, and weak 
poleward flow at depths that correspond to the near bottom portion of the water column 
on the shelf. The model results, however, do not seem to be as consistent with the ob- 
servations off Portugal, which have rather smaller mean meridional (v) velocities. 
b. Development of Eddies 

On day 50 the dynamic heights of the surface (relative to 2400 m) show a 
high pressure center (positive values over 2 cm) at approximately v= 640 kin, 60 km 
offshore (Figure 51 (left)). The high pressure region has a meridional extent that, with 


tine gets Jarger zonally and propagates offshore and to the south (compare Figure 5] 
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Figure 49. =‘ Vertical cross-shore section at day 50 of temperature (°C) at y=512 km 


for experiment 2 at day 50: Contour interval Is 1.0°C. 


(left) and (right)). Such propagation may be explained by the continuous offshore mo- 
tion in the upper layer, due to the uniform equatorward wind stress, and also due to the 
strong southward mean flow, especially on the eastern side of the ridge. Comparing the 
flow field of this experiment (day 50, after 30 days of equatorward wind) with experiment 
1 (day 20), the surface high pressure ridge developed at about the same alongshore po- 
sition, but in this case of uniform winds took a much longer time to generate eddies. The 
absence of negative wind stress curl explains that the poleward surface countercurrent 
offshore of the pressure ridge is much weaker and extends much less in the alongshiore 
direction. In experiment 1, the negative wind stress curl caused a much larger conver- 
gence and downwelling of isotherms on the west side of the pressure ridge, leading to 


much stronger poleward surface circulation offshore of the ridge, a much larger hor- 
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Figure 50. Vertical cross-shore section of alongshore-averaged meridional (¥) ve- 
locity (cm/s) for experiment 2 at day 50: Contour interval 1s 5.0 cms. 


Dashed lines denote equatorward velocities. 


zontal shear between the equatorward and poleward surface currents, and an earlier de- 
velopment of eddies due to a predominantly barotropic instability of the mean flow. 

By dav 75 (after 55 days of equatorward 3.5 dynes/cm? wind stress) the sur- 
face dynamic height field shows two anticyclonic circulations that develop in the region 
of the high pressure at positions of approximately y= 593 km and y=473 km, and one 
cyclonic circulation that develops within 50 km of the coast at v ~ 540 kin (Figure 51 
(right)). The surface isotherms for day 75 (Figure 52) indicate that the locations of both 
anucyclonic circulations correspond to the relatively warmer spots as seen by the 
concavities of the 16 and 15.5°C isotherms, and that the position of the cyclonic circu- 


lation corresponds to a region where colder isotherms extrude offshore. The cross- 
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Figure 51. Surface isopleths of dynamic height (cm) relative to 2400 m for exper- 
iment 2 at day 50 (left) and day 75 (right): Contour interval 1s 2.0 cm. 


Dashed lines denote negative values. 


stream derivatives of the time averaged (days 65-74) potential vorticity (Figure 53) show 
both zonal and vertical changes in sign in the upper ocean, near the eastern boundarv 
(within approximately 100 km of the coast), so, since the necessary conditions for 
barotropic and baroclinic instability are met, both mechanisms could be responsible for 


the generation of eddies in this experiment. 


At this point it is interesting to note that in experiment | (with negative 
wind stress curl) a much more intense eddy field with four anticyclonic eddies developed 
earlier (day 40) and no cyclonic eddies developed in the experiment. Such a response 
suggests that while in cases with no wind stress curl, the model may generate eddies with 
both cyclonic or anticyclonic rotations, for cases with forcing by winds with curl, the 
model is more likely to generate eddies with rotations consistent with the curl. 

From day 75 to day 90, the system of two anticyclonic and cyclonic circu- 
lations move southward. By day 90, from the meanders of the surface jet at y ~ 835 km, 
approximately 65 km offshore, a new anticyclone is developing (Figure 54 (left)). The 
pattern of this eddy field continues to intensify until the end of the experiment (day 105) 
and is displaced southwestward (Figure 54 (left) and (right)), probably due to advection 
by the mean flow. 

2. Experiment 3 (Forcing by a Time Series of the North/South Component of the 
Wind) 
a. General Picture of the Wind Time Series 

In experiment 3, the model was forced by a spatially uniform time series of 
the north/south component of the wind. The time series corresponds to the 6-hourly 
values of the meridional wind from | February to 31 December, 1986, computed from 
Fleet Numerical Oceanography Center synoptic surface atmospheric pressure analyses, 
interpolated to a three-degree latitude,longitude grid centered at position 39°N, 10°W. 

A summary of statistics related to the wind time series is shown in 
Table 4. The mean north/south wind for the series was -240 cm/s with a standard de- 
viation of 475 cms. The maximum poleward wind was 1490 cm/s and the maximum 
equatorward wind was -1430 cm/s. Monthly means of the north/south winds oscillated 
between 139 cm,s in December and -667 cm/s in July. February and December were the 
only months with a poleward mean wind. The stronger equatorward winds (mean values 
with magnitudes over 370 cm/s) occurred from May through August. September had a 
relatively small mean wind of -16 cm/s, but in October the mean of -338 cm/s indicated 
an increased influence of equatorward winds. By December, winds were poleward most 
of the tine and the mean value was 139 cm/s. The monthly standard deviations of the 
6-hourly north/south winds oscillated between a minimum of 256 cm's in October and 
a maximum of 497 cm's in February. The percentage of time with equatorward winds 
oscillated between a minimum of 39.2% in February and a maximum of 99.2% in July. 


Overall, this particular time series had prevailing and stronger equatorward mean winds 
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Surface isotherms (°C) for experiment 2 at day 75: Contour interval! 


1S ILO. 


Figure 52. 


(more than 75% of the time, and means larger than 370 cm/s) during the months from 


May through August. 
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Figure 53. Vertical cross-section of the cross-stream derivative of potential vorticity 
multiplied by the zonal grid size (°C / (im s)) for the time-averaged days 
63-74 of experiment 2: Values are scaled by 10°. Contour interval 1s 
0.1°C / (ms). Dashed contours denote negative values. The vertical 
cross-section was taken at v= 544 km. 
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iment 2 at day 90 (left) and day 105 (right): 


Dashed lines denote negative values. 


Table 4. MIONTHLY AND ENTIRE SERIES STATISTICS OF THE 
NORTH/SOUTH COMPONENT OF THE WIND VELOCITY 
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b. Response of the Flow to the Wind Forcing 
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Based on the structure of the time series, a few cases were selected and used 
to study the response of the flow. 
(1) Poleward IWinds 
(a) Fluctuating Poleward Winds During a Few Days in the 
Pre-upwelling Season— Here we focus on the time interval between day 5 (6 February) 
and dav 16 (17 February), in which the dominant poleward winds showed fluctuations 
with time scales of two or three days (Figure 55). Comparing the cross-sections of 
alongshore averaged meridional (v) velocitv for days 5 and 16 (Figures 56 and 57) we see 
that the existing equatorward surface current created by previous equatorward winds, 
eventually disappears and is replaced by a poleward surface coastal current, when the 
winds reverse. The zonal cross-sections of temperature (not shown) also demonstrate in 
a consistent fashion with the reversal of the flow, that the 15 and 16°C isotherms in the 


upper coastal ocean change from an upwelling to a downwelling situation. The response 
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of the model with a poleward current under the change of forcing to poleward winds 1s 
consistent with observations of mean currents over the central Oregon shelf during 
winter, When the wind is favorable for downwelling (Huver, 1983). 

(b) Short and Strong Pulses of Poleward Wind— A strong pulse 
of poleward wind between days 6! and 63 (3 - 5 April) with a peak velocity of about 10 
m's (after an existing situation of surface equatorward jet and poleward undercurrent), 
was used to evaluate the response of the flow. In this short period of time, the 
equatorward jet rapidly decreases in intensity and shrinks in cross-sectional area, the 
poleward undercurrent increases in intensity, expands in cross-sectional area, and dis- 
places the weak equatorward current offshore (see Figures 58 and 59). Notice that in 
comparing this case with the previous one, the equatorward current does persist, but is 
appreciably weakened, probably due to the less extensive time interval with poleward 
winds. 

(2) Equatorward Winds During the Upwelling Season | 

(a) Equatornwvard Winds with Intensification and Relaxation— A 
common feature of the prevailing equatorward winds during the upwelling season is the 
sequence of intensification and relaxation with time scales of 4-5 days. Intensification 
and relaxation are seen along the time record in intervals such as from day 64 to 78 (6 
to 20 April), from day 122 to 141 (3 to 22 June),from day 159 to 170 (10 to 21 July), and 
from dav 186 to 210 (6 to 30 August). 

The response of the system to intensification and relaxation 
may be exemplified in the interval between davs 135 and 141 (16 and 22 June). 
Figure 60 shows the time series of the meridional wind velocity for June, which includes 
the interval used in this example. On day 134 and 135 (16 June) the meridional wind 
velocities relax very close to zero. The wind then intensifies and reaches a maximum 
equatorward value of 12.1 m’s by the end of day 136 (17 June). Comparing the cross- 
sections of temperature at y= 624 km for day 135 (relaxation, Figure 61) and day 137 
(intensification, Figure 61) we notice that the upwelled 15 and 16°C isotherms respond 
near the surface with additional offshore displacement. which is consistent with increased 
offshore Ekman transport (Figure 63 (left) and (right)), as expected by the increased 
equatorward winds. The vertical structure of the alongshore averaged meridional (v) 
velocity also shows that, from day 135 up day 137, with the intensification of the wind, 
the surface jet increases in vertical extent and the cross-sectional area of the poleward 


undercurrent (values over 2.0 cms) 1s considerably reduced (compare Figures 64 and 65). 
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Figure 55. Time series of the 6-hourly values of the north/south component of the 
wind velocity (cm/s) for February 1986 at position 39°N, 10 °W:  Posi- 
tive and negative values represent poleward and equatorward winds, 


respectively. Day zero is | February 1986. 


The meridional (v) velocity in the core of the equatorward jet stays approximately the 
same. 

From day 137 (18 June) to day 141 (22 June), with the relaxa- 
tion of the wind, it is apparent that the equatorward surface current becomes shallower 
and the poleward undercurrent increases considerably in size and shoals (Figure 66). 
The 15 and 16°C isotherms recede toward the coast responding to the decreased offshore 
Ekman transport, due to the relaxation (compare Figure 67 with Figure 62), but the 
configuration of the isotherms in the upper 200 m stays upwelled. Off Oregon, under 
variable winds during the upwelling season, observations have shown that isopycnals 
remain tilted upward toward the coast while the near-surface density field changes rap- 


idly with the wind (Huver, 1988). Observed motion of light (warmer) surface waters 
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Figure 56. Vertical cross-shore section of alongshore-averaged meridional (v) veloc- 


ity (cm/s) for experiment 3 at day 5: Contour interval is 2.0 cms. 


Dashed lines denote equatorward velocities. 


onshore when equatorward winds relax and offshore when winds intensify ts consistent 
with the model response just described. An experiment conducted during June and July 
of 1976, off Portugal in the region of the Nazaré canyon (latitude 39.5°N), using an 
airborne radiation thermometer to survev sea surface temperature, showed that with the 
slackening of the northerly winds, the thermal pattern changed consistently with a 
shoreward penetration of Warmer oceanic water (Fiuza, 1982). A sequence of thermal 
infrared satellite rmages between 31 August and 2 September 1981, under the action of 
upwelling favorable winds, showed a continuous band of upwelled waters at the north 
coast of Portugal, as far south as Cape Roca (latitude 38.6°N), an isolated focus of up- 


welling on the southwest coast, and offshore (over the slope) a meridional band with 
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Figure 57. Vertical cross-shore section of alongshore-averaged meridional (v) veloc- 
ity (cm/s) for experiment 3 at day 16: Contour interval is 2.0 cm/s. 


Dashed lines denote equatorward velocities. 


large zonal gradients of temperature. During the subsequent relaxation of the 
equatorward winds, the satellite data for 3,4, and 10 September 1981 indicated the 
gradual disappearance of the upwelling focus, the weakening of the thermal gradients, 
and a progressive shoreward advance of warmer oceanic waters (Fiuza, 1984). Thus, the 
results that we get in experiment 3, under relaxation of the upwelling favorable winds, 
seem consistent with observations off both Oregon and Portugal. 

(b) Quasi-steady Equatorward Winds During a Few Days After 
Short Periods of Poleward Winds— A situation that (in this particular year) is repeated 
a few times, consists of moderate quasi-steady equatorward winds of approximately 5 


m's that persist for about 6-9 days, preceded by weak poleward winds. Examples may 
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Figure 58. Vertical cross-shore section of alongshore-averaged meridional (¥) veloc- 
ity (cm/s) for experiment 3 at day 61: Contour interval is 2.0 cm,s. 


Dashed lines denote equatorward velocities. 


be seen in the time series (Figure 3) in the intervals from day 30 to 37 (3 to 10 March), 
from day 45 to 52 (18 to 25 March), from day 55 to 60 (28 March to 2 April), from day 
82 to 91 (24 April to 3 May), from dav 97 to 102 (8 to 14 May) and from day 256 to 264 
(15 to 23 October). Such situations in this particular series seem to be more frequent in 
the early upwelling season (March and April). The response of the flow is very similar 
to that described previously for the intensification. That is, the equatorward surface 
current near the coast tends to expand in cross-section and the poleward undercurrent 
is reduced or even disappears, as shown in Figure 68 and Figure 69 which represent the 
beginning (dav 96 or 8 Mav) and end (day 102 or 14 Mav) of the interval with quasi- 


steady, moderate, equatorward winds. 
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Figure 59. Vertical cross-shore section of alongshore-averaged meridional (¥) veloc- 
ity (cm/s) for experiment 3 at day 63: Contour interval is 2.0 cm/s. 


Dashed lines denote equatorward velocities. 


In summary, with increased equatorward winds, the surface 
equatorward jet responds by expanding its cross-sectional area and increasing its velocity 
while the poleward undercurrent responds by reducing its cross-sectional size and de- 
creasing its velocity. Relaxation of equatorward winds causes the poleward undercurrent 
to expand and shoal. Velocities measured on a zonal section off Point Sur, California 
in November 1988 showed that the inner shelf poleward flow along the coast can be 
weak to nonexistent during strong upwelling, but strengthens and propagates seaward 
during relaxations (Reece, 1989). In a consistent fashion, the model responds to inten- 
sified upwelling favorable winds by submerging the poleward undercurrent, and responds 


to wind relaxations by allowing the poleward current to shoal. 


81 


1000.0 
Sidtige Gl 26 aa 


~1000.0 


en, 
Za) 
os 
E 
hy 
~ 
— 
O 
2 
® 
> 
22 
= 
= 
ns 
ors 
— & 
(‘= 
ey 
(S 
Oo 
Q. 
e 
OQ 
he) 
® 
asf 
e) 
es 
Ww) 
oO 
Cc 
2 
< 


—2000.0 
hi-}------------5----' 





Figure 60. 


eels Rel 


ieee lle tts Celt ee 


ma mm eet ee nh - - - - - - - ee ee ye ee ee = 
Se ee 


let ole 


i, co 


Gepe == see. 


14 | 18 20 
Time in days 


Time series of the 6-hourly values of the north/south component of the 
Wind velocity (cm/s) for June 1986 at position 39°N, 10 °W: Positive 
and negative values represent poleward and equatorward winds, re- 


spectively. Day zero is | June 1986 (day 120 in experiment 3). 


(3) Development of Eddies During the Upwelling Season. On day 173 (24 


July), as seen in the held of dynamic height of the surface relative to 2400 m (Figure 70 


(right)), there is a meandering of the surface jet between positions y ~ 640 km and y ~ 


896 km, which relates and is in phase with the thermal surface signal shown for the same 


date in Figure 70 (left). The estimated wave length for the meander 1s 125 km. 


By day 188 (8 August) the meander has increased in amplitude, a 


closed anticyclonic circulation develops at y ~ 748 km, and the surface temperature 


pattern shows extended offshore cold water filaments corresponding to the cyclonic cir- 


culations of the meandering jet at y ~ 720 km and y ~ 841 km ( Figure 7I(left) and 


(right)). The cross-stream derivatives of the time averaged (days 180-189) potential 


vorticity (Figure 72) show a change in sign from positive above to negative below, 


T FIELD at Day 135 at y = 624 km 
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Figure 61. Vertical cross-shore section of temperature (°C) at y= 624 kim for ex- 


periment 3 at day 135: Contour interval is 1.0°C. 


across a horizontal line at a depth of about 35 m, which corresponds to a necessary 
condition for baroclinic instability. Therefore a possible mechanism for the generation 
of the eddies during the upwelling season is the baroclinic instability due to the vertical 
shear associated with the equatorward surface jet and poleward undercurrent. 

On day 197 (17 August), as seen in the field of dynamic heights of the 
surface relative to 2400 m (Figure 73 (right)) the anticyclone has expanded in size (as 
compared with day 188) and two new closed cyclonic circulations show up at y ~ 832 
km and y ~ 704 km. The closed anticyclonic circulation is also identified by the zonal 
cross-section at y= 736 km representing the vertical structure of the temperature and 
meridional (v) velocity on dav 197 (Figures 74 and 75). In the position of the warm core 


eddy (centered at about 64 km offshore) the 17°C isotherm is displaced from the surface 
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Figure 62. Vertical cross-shore section of temperature (°C) at y= 624 km for ex- 


periment 3 at day 137: Contour interval is 1.0°C. 


downward as much as 38 m and rises back to the surface. Consistent with this thermal 
structure, the meridional (v) velocity is towards the north (larger than 2 cm s in the up- 
per 100 m) on the west side of the eddy and southward (larger than 2 cm/s in upper 300 
in) on the east side, where it interacts with the southward coastal jet. The northern and 
southern closed cyclonic circulations can also be recognized by looking at the structure 
of the isotherms and meridional (v) velocities in the vertical zonal sections across the 
eddies. For the southern cvclone, at v= 704 km, the configuration of the coastal 
isotherms in the upper 300 m shows a cold dome (raised isotherms) in the central posi- 
tion of the eddv (Figure 76). As seen in Figure 77, the distribution of the meridional 
(v) velocity is indicative of a southward current on the west side of the cyclone (between 


24 and 74 km from the coast) and a northward current on the east side of the cyclone 
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Figure 63. Surface isopleths of zonal (u) velocity (cm/s) for experiment 3 at day 135 
(left) and day 137 (right): Contour interval is 5.0 cm/s. Dashed lines 


denote offshore velocities. 


(within 24 km from the coast). The vertical extent of both cyclonic eddies 1s estimated 
to be about 300 m. In between a dipole pair of eddies composed of the anticyclone and 
the northern cyclone (Figure 73 (right)), lies a cold water filament identified by the rel- 
ative maximum offshore displacement of the 16.5°C isotherm, that at y ~ 705 km gets 


as far as 62 km from the coast (Figure 73 (left)). 
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Figure 64. 
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Vertical cross-shore section of alongshore-averaged meridional (¥) veloc- 
ity (cm/s) for experiment 3 at day 135: Contour interval is 2.0 cms. 


Dashed lines denote equatorward velocities. 


By day 213 (2 September), a second closed anticyclonic circulation 1s 


apparent in the field of dynamic height of the surface relative to 2400 m at approxi- 


mately v= 848 km (Figure 78 (left)). The structures of the temperature and meridional 


(v) velocities in the zonal vertical sections (not shown) are also indicative of a warm core 


ring with an antucyclonic circulation centered at about 77 km offshore. The vertical ex- 


tent of the anticyclonic eddy, as seen from the meridional (v) velocity and temperature 


structures (not shown), is about 450 m. 


From dav 213 to day 225 (from 2 to 14 September), each of the 


cyclonic circulations splits in two and results in four closed cyclonic circulations, that 


span the coast with the donunant axis along the coast (Figure 78 (left) and (right)). The 
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Figure 65. Vertical cross-shore section of alongshore-averaged meridional (v) veloc- 
ity (cm/s) for experiment 3 at day 137: Contour interval is 2.0 cmys. 


Dashed lines denote equatorward velocities. 


cross-section of the alongshore averaged meridional (v) velocity (not shown) shows a 
poleward nearsnore surface current within about 30 km of the coast, with speeds above 
2 cms in the upper 200 m. Offshore of 30 km, a southward surface current is evident. 
From day 225 to day 270 (from 14 September to 29 October), the two anticyclones 
propagate southwestward, probably due to a combination of advection by the mean flow 
and #-induced motion. 

(4) End of the Upwelling Season. The cyclonic circulations that spanned 
the coast become almost imperceptible by dav 270 (Figure 79 (left)). The coastal cir- 
culation consists of a weak equatorward current offshore and a weak poleward surface 


current nearshore. The weakening of the eddy field near the coast may correspond to a 
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Figure 66. Vertical cross-shore section of alongshore-averaged meridional (¥) veloc- 
ity (cm/s) for experiment 3 at day 141: Contour interval 1s 2.0 cm/s. 


Dashed lines denote equatorward velocities. 


real decrease of energy in the mesoscale features by the end of the upwelling season. 
During November and December, the alongshore winds, with decreased mean 
equatorward velocities and periods of moderate or strong poleward velocities (9 to 10 
November, 12 to 13 November, 27 November to 7 December), have decreased the 
equatorward surface current and probably weakened the eddies by modifying the ther- 
mal structure. By the end of the experiment on day 332 (30 December), the two 
anticyclonic eddies no longer exist, three very weak cyclonic circulations are still evident, 
and the surface equatorward current is now much less intense and has a smaller cross- 


sectional area than during the peak of the upwelling season (Figure 79 (right)). 
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Figure 67. Vertical cross-shore section of temperature (°C) at y= 624 km for ex- 


periment 3 at day 141: Contour interval is 1.0°C. 
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Figure 68. Vertical cross-shore section of alongshore-averaged meridional (¥) veloc- 
ity (cm/s) for experiment 3 at day 96: Contour interval is 2.0 cm/s. 


_ Dashed lines denote equatorward velocities. 
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Figure 69. Vertical cross-shore section of alongshore-averaged meridional (v) veloc- 


ity (cm/s) for experiment 3 at day 102: Contour interval is 2.0 cm/s. 


Dashed lines denote equatorward velocities. 
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Figure 70. Surface isopleths of (left) temperature (°C) and (right) dynamic height 


(cm) relative to 2400 m for experiment 3 at day 173: Contour interval 


is 0.5°C for temperature and 1.0 cm for dynamic height. Dashed lines 


denote negative values. 
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Surface isopleths of (left) temperature (°C) and (right) dynamic height 


Figure 71. 
Contour interval 


(cn) relative to 2400 m for experiment 3 at day 188: 
1s 0.5°C for temperature and 1.0 cm for dynamic height. Dashed lines 


denote negative values. 
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Figure 72. Vertical cross-section of the cross-stream derivative of potential vorticity 
multiplied by the zonal grid size (°C / (m s)) for the time-averaged days 
180-189 of experiment 3: Values are scaled by 10°. Contour interval 
is 0.1°C / (ms). Dashed contours denote negative values. The vertical 
cross-section was taken at v= 752 km. 
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Surface isopleths of (left) temperature (°C) and (right) dynamic height 


Figure 73. 
Contour interval 


(cm) relative to 2400 m for experiment 3 at day 197: 
is 0.5°C for temperature and 1.0 cm for dynamic height. Dashed lines 
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denote negative values. 
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Figure 74. Vertical cross-shore section of temperature (°C) at y = 736 km for ex- 


periment 3 at day 197: Contour interval is 1.0°C. 
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Figure 75. Vertical cross-shore section of meridional (v) velocity (cm/s) at y = 736 


km for experiment 3 at day 197: Contour interval is 2.0 cm/s. Dashed 


lines denote equatorward velocities. 
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Figure 76. Vertical cross-shore section of temperature (°C) at y = 704 km for ex- 


periment 3 at day 197: Contour interval is 1.0°C. 
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Figure 77. Vertical cross-shore section of meridional (v) velocity (cm/s) at y = 704 


kin for experiment 3 at day 197: Contour interval is 2.0 cm,s. Dashed 


lines denote equatorward velocities. 
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Figure 79. Surface isopleths of dynamic height (cm) relative to 2400 m for exper- 
iment 3 at day 270 (left) and day 332 (right): Contour interval is 1.0 


cm. Dashed lines denote negative values. 
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VY. SUMMARY AND RECOMMENDATIONS 


This study used a high-resolution, multi-level, primitive equation ocean model to 
investigate the role(s) of climatological and time dependent wind forcing as possible 
generation mechanisni(s) for the coastal jets, undercurrents, fronts, and eddies in the 
eastern boundary current region off the west coast of the Iberian Peninsula. A band of 
steady equatorward winds, uniform alongshore but with zonal variability (a large nega- 
tive wind stress near the coast, and narrow band of negative wind stress curl followed 
by regions of positive and negative wind stress curl offshore) was used as forcing in ex- 
periment 1. A band of spatially uniform, poleward winds fromm day 0 to day 20, followed 
by a reversal to spatially uniform, equatorward winds from days 21 to 105 was used as 
forcing in experiment 2. A band of spatially uniform, time dependent winds corre- 
sponding to the north/south component of the 6-hourly winds, computed from synoptic 
surface pressure analyses for the period fron. 1 February to 31 December 1986, was used 
as forcing in experiment 3. All experiments were performed on a f-plane, in an idealized, 
flat-bottom oceanic regime along the eastern boundary off Portugal. 

On day 5 of experiment 1, in the meridional band of negative wind stress curl, there 
was a change in sign of the zonal gradients : of surface zonal velocity, of surface dynamic 
height, and of the upper ocean temperature field. The upwelled isotherms near the coast 
Were consistent with vertical motion near the coast to compensate the Ekman layer d1- 
vergence due to the equatorward winds. Further offshore, in the region of negative wind 
stress curl (due.to the decreased equatorward winds), the westward Ekman transport was 
reduced, resulting in surface laver convergence and downwelling. As expected, the dis- 
tribution of the surface dynamic heights relative to 2400 m (which are proportional to 
temperature) shows a meridionally extended ridge in the band of negative wind stress 
curl. The equatorward surface current nearshore and poleward surface current offshore 
of the surface pressure ridge are approximately in geostrophic balance with the pressure 
gradient forces due to the horizontal distribution of density. The flow structure on day 
5 was consistent with a similar experiment conducted by McClain ez al. (1986) for the 
north and west coast of Galicia, in which a maximum southward speed resulted 
shoreward of the wind shear zone and a maximum northward speed occurred offshore 
of the zone. In our experiment 1, the maximum equatorward velocity was larger by a 


factor of two, which may be explained by the use of a larger equatorward wind stress 


at the coast (3.5 dynes/cn? versus 3.0 dynes/cn? ) and a larger negative wind stress curl 
offshore. 

By day 20 of experiment |, both equatorward and poleward surface currents were 
stronger, and a poleward undercurrent was present adjacent to the boundary at depths 
between ~200 and 500 m, consistent with previous results of McCreary et al. (1987) and 
Batteen et al. (1989). There was evidence of eddies at y ~735 km and y ~191 km super- 
imposed on this system. 

By day 40 of experiment 1, four, well developed anticyclonic eddies were present 
along a meridional axis located close to the band of negative wind stress curl. All eddies 
were likely generated due to the anticyclonic shear of the wind in this region and each 
had a thermal structure consistent with warm core rings. Four cold water filaments were 
lying in places where the circulation of each anticyclone advected the upwelled water 
offshore. The alongshore wave length of the eddies was estimated, based on spectral 
analysis, to be approximately 154 km. 

By day 60 of experiment 1, the two northern and the two southern eddies had 
merged/coalesced into one northern and one southern eddy. From days 60 to 90, the two 
eddies propagated westward and expanded in size. The estimated westward propagation 
rates Were 1.4 km, day and 2.3 km,day for the northern and southern anticyclones, re- 
spectively. The faster rate for the southern eddy is qualitatively consistent with a faster 
B-induced motion at lower latitudes. 

By dav 90 of experiment I, the eddies had stabilized in size and, as estimated bv the 
spectral analysis, had a peak energy at a dominant wave length of ~300 km. From day 
90 to the end of experiment | (day 154), the two anticvclones continued to propagate 
westward. By day 150, two new anticyclones were generated in the same region of neg- 
ative wind stress curl as before. The thermal and velocity structures were consistent via 
the thermal wind relation. 

The equatorward nearshore surface current and poleward surface countercurrent 
offshore, reproduced early in the experiment, Were qualitatively consistent with a similar 
experiment (McClain et a/., 1986) and with northward transport near the shelf break, as 
indicated by Coastal Zone Color Scanner (CZCS) imagery off the west coast of Galicia. 
Hlowever, due to the lack of available observations, it is not known if the instability of 
the flow and the resulting mesoscale eddies in this experiment correspond to reality off 
the IP. While the negative wind stress curl has been observed off the west coast of 
Galicia from 18 to 23 April 1982, and was probably orographically induced by the rias, 


no observations exist to verify how far downstream of the Galicia coast this negative 
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wind stress curl may occur. With the lack of observations, the results of experiment | 
may best be considered a process oriented study designed to study the influence of a 
small-scale, large negative wind stress curl superimposed on a large scale climatological 
positive wind stress curl offshore in an eastern boundary current sysyvem. 

Experiment 2 was conceived to study a “spring transition”. On day 20 (end of the 
period of 1.2 dynes/cm? poleward wind stress) the flow responded with onshore Ekman 
transport convergence in the Ekman layer, and downwelling near the coast. The surface 
dynamic heights relative to 2400 m consistently indicated increasing values towards the 
coast and, as expected, a poleward surface current extended along the eastern boundary 
with hghter/warmer waters at the coast. Observations over the shelf off Oregon during 
the winter (prevailing poleward winds) have shown poleward flow over the shelf without 
mean vertical shear (Huyer, 1983). In experiment 2 the poleward surface current 
(alongshore averaged maximum meridional velocity of ~ 12 cm/s extending down to 
~150 m) showed vertical shear, which may be explained via the thermal wind relation 
by the exaggerated downward displacement of the isotherms in the coastal upper ocean. 
There was also a Weak equatorward undercurrent. 

Experiment 2 showed a relatively fast (1-2 davs) response of the svstem to a reversal 
to equatorward winds of 3.5 dynes/cn? . On day 22 (2 days after the transition to 
equatorward wind stress), the model responded with offshore Ekman transport tn the 
upper laver, upwelled isotherms near the coast, and the replacement of the previous 
poleward surface current by an equatorward surface current. On day 30, with increasing 
time under the influence of the equatorward winds, the upwelled 16 and [5°C isotherms 
intercepted the surface and were displaced offshore and the surface equatorward coastal 
current, associated with the front in the upper ~70 mm, had a jet-like structure which 
showed large vertical shear in the alongshore velocity. The poleward undercurrent, with 
an averaged alongshore velocity of ~5 cm/s, was present nearshore between ~320 and 
450 m depth. In experiment 2, the changes of the thermal and velocity structures near 
the coast were consistent with observations of the 25 through 31 March 1975 spring 
transition off Oregon (Iluyer, 1983) by having upward sloping isotherms, a mean 
southward surface current, and strong vertical shear. No observations are available to 
document a “spring transition” off Portugal. 

From day 30 to day 50 the surface isotherins were continuously displaced offshore, 
the surface equatorward jet widened as far as 64 km from the coast (day 50) and became 
shallower inshore. [The upper portion of the poleward undercurrent became closer to the 


surface (depth of ~120 m at the coast) and a weak surface poleward countercurrent de- 
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veloped offshore of the equatorward jet in the middle model domain. The southward 
surface jet and poleward undercurrent at depths below 120 m were consistent with the 
typical summer coastal upwelling regime off Oregon, as described by Huver (1983). OfF 
Portugal, current meter observations (Fiuza, 1984) from 12 June to 7 July 1980 (peak 
of the upwelling season) on the inner and outer shelf at latitude 38° 06’N showed weak 
mean southward flows (-1.7 and -2.1 cm/s for the inshore and near-shelf break stations, 
respectively). Hlowever, these observations were taken under a regime of quasi- 
permanent southward, but fluctuating winds, and topographic features, such as the 
Setubal embayment and Setubal canyon in the vicinity of the current meter stations, 
thus making comparisons with the model results difficult, because the model assumes a 
Straight coastline, flat bottom and was forced by a steady, relativelv large wind stress. 

Two anticyclonic eddies were generated between days 70 and 75 of experiment 2 at 
positions y ~593 km and y ~473 kmina region of high surface dynamic height that was 
existent in approximately the middle of the model domain, 60 km offshore. One closed 
cyclonic circulation also existed on dav 75, that developed from the meandering jet at v 
~ 540 km near the coast. From day 75 to the end of the experiment (day 105), a new 
anticvclone developed from meanders of the coastal jet and was present on day 90 at v 
~835 km, approximately 65 km offshore. The eddies intensified and were displaced con- 
tinuously southwestward probably due to advection by the mean flow. Experiment 2, 
like experiment 1, can be considered a process-oriented study to give insight on the 
coastal ocean processes associated with “spring transitions” along eastern boundaries. 

The time series of the wind velocities used to force the model in experiment 3 sug- 
gests that pulses of strong or moderate poleward winds with ume scales of a few days 
(typically ~ 4-6 days) may be important, especially during the pre-upwelling season 
(February) and also after the decay of upwelling (November and Deceinber). In exper- 
iment 3, between days 5 and 16 (6 to 17 February), the existing equatorward coastal 
surface current, created by previous equatorward winds, disappeared and was replaced 
by a poleward surface coastal current. Situations of short pulses of strong poleward 
winds (for instance from day 61 to 63) may not be enough to completely reverse a pre- 
existing equatorward surface current, but the current tends to decrease in intensity and 
is reduced in cross-sectional area, while the poleward undercurrent increases in intensity, 
expands in cross-sectional area, and displaces the weakened equatorward current off- 
shore. 

During the upwelling season, fluctuations of the prevailing equatorward winds with 


time scales of a few days (typically 4-5 days) correspond to events of intensification and 
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relaxation. With intensification of the winds favorable for upwelling, the upwelled 
isotherms near the surface respond with additional offshore displacement consistent with 
increased offshore Ekman transport, the surface equatorward jet increases in vertical 
extent, and the poleward undercurrent 1s considerably reduced. With the relaxation of 
the equatorward winds, the surface isotherms tend to recede towards the coast, but the 
configuration of the isotherms in the upper 200 m stays upwelled. Off the coastal of 
Portugal, surveys of sea surface temperature using an airborne radiation thermometer 
(Fiuza, 1982) and satellite infrared imagery (Fiuza, 1984) have indicated onshore 
progression of warmer oceanic waters in association With relaxation of the northerly 
winds. During the relaxation, the equatorward surface current became shallower and the 
poleward undercurrent got closer to the surface and increased in size. Experiment 3 also 
showed that by the middle and end of the upwelling season, mesoscale eddies are likely 
to be generated off the west coast of Portugal. On day 188 (8 August) from a meander 
of the surface coastal jet, a closed anticyclonic circulation developed at y ~ 748 km. 
On day 197 (17 August) two new closed cyclonic circulations show up at v ~ 832 km 
and y~ 704 km. The thermal structure of the anticvclone shows downward displacement 
of the upper ocean isotherms in the center of the eddy. The thermal structures of the 
cyclonic eddies show raised upper ocean isotherms in the center of the eddies. Distrib- 
uuons of the velocities associated with the eddies are consistent with the thermal struc- 
tures via the thermal wind relation. A cold water filament is recognized at the surface 
on dav 197 (17 August), Iving between a dipole pair of eddies composed of the 
anticyclone and the northern cyclone. On day 213 of experiment 3 (2 September), a 
second anticyclone was present at ~848 km. The two cyclonic eddies positioned closer 
to the coast than the anticyclone split, resulting in four cyclones that spanned the coastal 
region by day 225 (14 September). During November and December, with the decreased 
influence of the winds favorable for upwelling and with some time intervals of 
strong’moderate poleward winds (days 281 to 283, 284 to 286, and 299 to 310), the 
eddies, as recognized in the distribution of surface dynamic height, weakened and dis- 
appeared. On day 332 (30 December) the two anticyclones were no longer evident and 
only weak remains of the cyclonic circulations were superimposed on the also weakened 
surface coastal current svstem. 

The results from the experiments support the hypothesis that wind forcing, and 
particularly time dependent winds, can be a significant generation mechanism for eddies, 
jets, fronts, and upwelling filaments off the west coast of the IP. It should be noted, 


however, that this study employed the constraints of a regular, straight coastline and a 
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flat bottom to isolate and examine the effects of wind forcing. Future studies should use 
climatological, spatially resolved winds (values on a grid with points one-degree or even 
less apart) for the region and also time dependent, spatially resolved winds. Features like 
the irregular coastline and bottom topography should also be included in the numerical 
Studies, because they are likely to be important mechanisms in controlling the locations 
of mesoscale features in the study region. Inclusion of salinity, which is important in 
controlling the density at the depths of the Mediterranean Water influence, may also 
improve the accuracy of the results in the experiments. Besides the numerical studies, 
observational programs on the shelf and slope should be implemented including current 
meter arrays, Lagrangian current drifters, remote sensing, and CTD surveys with ap- 
propriate space/time resolution to synoptically map mesoscale structures. Such a pro- 


gram would allow realistic comparisons with the numerical studies. 
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